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Abstract 

The Standard Model Higgs sector, extended by one weak gauge triplet of scalar 
fields with a very small vacuum expectation value, is a very promising setting to 
account for neutrino masses through the so-called type II seesaw mechanism. In this 
paper we consider the general renormalizable doublet /triplet Higgs potential of this 
model. We perform a detailed study of its main dynamical features that depend on 
five dimensionless couplings and two mass parameters after spontaneous symmetry 
breaking, and highlight the implications for the Higgs phenomenology. In particular, 
we determine i) the complete set of tree-level unitarity constraints on the couplings 
of the potential and ii) the exact tree-level all directions boundedness from below 
constraints on these couplings. When combined, these constraints delineate precisely 
the theoretically allowed parameter space domain within our perturbativc approxi- 
mation. Among the seven physical Higgs states of this model, the mass of the lighter 
(heavier) CVeven state h^ {H^) will always satisfy a theoretical upper (lower) bound 
that is reached for a critical value //c of fi (the mass parameter controlling triple 
couplings among the doublet/triplet Higgses). Saturating the unitarity bounds we 
find an upper bound m/jO < 0(0.7 — ITeV), while the upper bound for the remaining 
Higgses lies in the several tens of TeV. However, the actual masses can be much 
lighter. Wc identify two regimes corresponding to jU > /i^ and fi < fic- In the first 
regime the Higgs sector is typically very heavy and only h^ that becomes SM-like 
could be accessible to the LHC. In contrast, in the second regime, somewhat over- 
looked in the literature, most of the Higgs sector is light. In particular, the heaviest 
state becomes SM-likc, the lighter states being the CVodd Higgs, the (doubly) 
charged Higgses and a decoupled possibly leading to a distinctive phenomenology 
at the colliders. 

* corresponding author 



1 Introduction 



One of the major goals of the LHC is to uncover the mechanism underlying the electroweak 
symmetry breaking and thereby the origin of the weak gauge boson and fermion masses. 
Moreover, observation of neutrino oscillations has shown that neutrinos are massive (for a 
review see for instance [1] and references therein). Such masses do not necessarily require 
physics beyond the standard model (SM), since one can accommodate a (Dirac) mass 
through a Yukawa coupling assuming a right-handed neutrino similarly to the other massive 
fermions. However, the introduction of such a right-handed state, whose only role is to allow 
for non-zero neutrino masses while being neutral under all the SM interactions, might seem 
rather mysterious. Furthermore, in contrast with the other right-handed fermion states 
of the SM, the right-handed neutrino allows also for a Majorana mass that is invariant 
under the SM gauge group but violates lepton number. These features make plausible 
the existence of new flavor physics beyond the SM associated with the neutrino sector. 
Probably, one of the most attractive aspects is the ability to induce naturally the tiny 
neutrino masses from this new flavor physics sector [2] . The celebrated seesaw mechanism 
[31 m [5] relating directly the smallness of the neutrino masses to the presence of a large 
new scale A through m^, ~ f^/A, when A ^ v where v denotes the electroweak scale, is 
realized in a grand unified context (GUT) comprising right-handed neutrinos and often 
dubbed type I seesaw. It can also be achieved without right-handed neutrinos through an 
extended Higgs sector including an SU{2)l triplet scalar field, type II seesaw [61171181191 flO] . 
or by including two extra matter multiplets in the adjoint of SU{2)l, type III seesaw [llj . 
or a hybrid type mixture of type I and type III [12], |T3l HH [15] . 

If such extended sectors are too heavy to be directly accessible to TeV scale experiments, 
they could still be indirectly probed through distinctive low energy effective operators in 
the neutrino sector [TB]. In the present paper we will rather focus on the possibility 
of accessing directly the Higgs sector per se of the type II scenario, studying general 
dynamical constraints which originate from the potential that couples the Higgs doublet 
and the Higgs triplet. However, given the present theoretical uncertainties, we do not 
commit to any specific GUT or flavor physics scenarios beyond the SM. In particular 
mass parameters such as /i and Ma will not necessarily take large GUT scale values, even 
though such a configuration is included in the analysis. We will even consider regimes 
with very small fi (^ Gp ). As noted in [17], such a small /i makes all the Higgs sector 
accessible to the LHC. Here we carry out a complete study, taking into account the full 
set of renormalizable operators present in the potential. The aim is to exhibit the various 
possible regimes consistent with the dynamical constraints dictated by the potential and 
their consequences on the phenomenology of the extended Higgs sector. Most of these 
operators are often neglected in the existing phenomenological studies of type II seesaw, 
based on the fact that after spontaneous symmetry breaking their effects are suppressed 
by the small Higgs triplet vacuum expectation value (VEV), Vt, when compared to the 
electroweak scale. This is, however, not justified when studying the small n regimes just 
mentioned, where fi can be of order vt- In this case the detailed dynamics leads to an 
interesting structure of the Higgs sector. 

The paper is organized as follows: in section 2, we present the ingredients of the model, 
the physical Higgs states and mass spectrum, as well as a parameterization of the po- 
tential parameters in terms of the physical masses. In section 3, we discuss some of the 
phenomenological and theoretical constraints on the parameters related to precision mea- 
surements, the absence of tachyonic Higgs modes, as well as the presence of false vacua. In 
section 4, we provide a thorough study of the boundedness from below of the potential and 
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establish for the first time simple necessary and sufficient conditions on the couplings that 
are valid for all field directions. The unitarity constraints are analyzed in detail in section 
5, through the study of all the scalar scattering channels. In section 6, we combine in an 
analytical compact form the constraints obtained in sections 4 and 5. Section 7 presents the 
behavior of the CVeven Higgs masses as functions of the potential parameters, highlighting 
theoretical upper and lower mass bounds and identifying different regimes that give better 
insight into the overall Higgs sector phenomenology, as well as the determination of unitar- 
ity mass bounds on the lightest Higgs. Section 8 is devoted to a short review of the salient 
features of the Higgs phenomenology at the colliders as well as to specific illustrations of 
our results. We conclude in section 9 and give some technical details in the appendices. 



2 The model 

We start by recalling the scalar potential and the main properties of the Higgs physical 
eigenstates after EWSB as well as the corresponding eigenmasses and mixing angles. We 



give the expressions without neglecting any of the couplings appearing in Eq. (2.4) nor 
making any specific assumption about the magnitudes of fi, m^j and M\ which would 
originate from the unknown underlying high energy theory. The results of this section fix 
the notations and will serve for the completely model-independent analysis carried out in 
the subsequent sections. 



2.1 The Higgs potential 

The scalar sector consists of the standard Higgs weak doublet H and a colorless scalar field 
A transforming as a triplet under the SU{2)l gauge group with hypercharge Ia = 2, so 
that H ~ (1, 2, 1) and A ~ (1, 3, 2) under the 5[/(3)c x SU{2)l x f/(l)y. 

Under a general gauge transformation U{x), H and A transform as if — )■ U{x)H and 
A — 7- U{x)AU'^{x). One can then write the most general renormalizable and gauge invariant 
Lagrangian of this scalar sector as follows: 

C = {D,H)\D^'H) + Tr{D,A)\D^A) - V{H, A) + ^Yukawa (2.1) 
where the covariant derivatives are defined by 

D,H = d,H + igT'^W'^^H + i^^B^H (2.2) 

= d,A + ig[T'^W;, A] + ig'^B.A (2.3) 

(VF^, g), and (5^, g') denoting respectively the SU{2)l and U{1)y gauge fields and cou- 
plings and = cr"/2, with a" (a = 1,2,3) the Pauli matrices. The potential ^(ii. A) is 
given by, 

V{H, A) = -mlH^H + ^{H^ Hf + MlTr{A'^ A) + [fi{H^ia''A^ H) + h.c] 

+Xi{H^H)Tr{A'^A) + XiiTrA'^A)^ + X^TriA'^AY + X^H^ AA^ H (2.4) 

where Tr is the trace over 2x2 matrices. /^Yukawa contains all the Yukawa sector of the SM 
plus one extra Yukawa term that leads after spontaneous symmetry breaking to (Majorana) 
mass terms for the neutrinos, without requiring right-handed neutrino states. 
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vCvukawa 3 -Y^L'^C ® ia^AL + h.c. (2.5) 

where L denotes SU{2)l doublets of left-handed leptons, denotes neutrino Yukawa 
couplings, C the charge conjugation operator, and we have suppressed falvor indices for 
simplicity. Although part of the type II seesaw model, we will refer to the above model 



Eq. (2.1) as the doublet-triplet-Higgs-Model (DTHM) since in this paper we are mainly 
interested in the scalar sector, bringing up only occasionally the content of the Yukawa 
sector ^Yukawa and the related neutrino masses issue. 

Defining the electric charge as usual, Q = h + ^ where / denotes the isospin, we write 
the two Higgs multiplets in components as, 

where we have used for convenience the 2x2 traceless matrix representation for the triplet 



The potential defined in Eq. (2.4) exhausts all possible gauge invariant renormalizable 
operators. For instance a term of the form Asif^A^Aif , which would be legitimate to add 
if A contained a singlet component, can actually be projected on the Ai and A4 operators 



appearing in Eq. (Q thanks to the identity WA'^AH + WAA^H = WHTr{A'^A) which 
is valid because A is a traceless 2x2 matrix. This simply amounts to redefining Ai and 
A4 such as Ai + A5 — 7- Ai, A4 — A5 — > A4. The potential thus depends on five independent 
dimensionless couplings A and Xi,{i = 1,...4) and three mass parameters, m'jj,M^ and 
fi. In the present paper we will assume all these parameters to be real valued. Indeed, 
apart from the fi term, all other operators in V are self-conjugate so that, by hermiticity 
of the potential, only the real parts of the A's and the m^, mass parameters will be 
relevant. As for /i, the only parameter that can pick up a would-be CP phase, this phase is 
unphysical and can always be absorbed in a redefinition of the fields H and A. One thus 
concludes that the DTHM Lagrangian is CP conserving (see also the discussion in [18]). 
Moreover, V depends on five complex (or ten real) scalar fields. 

Assuming that spontaneous electroweak symmetry breaking (EWSB) is taking place at 
some electrically neutral point in the field space, and denoting the corresponding VEVs by 



one finds after minimization of the potential Eq.(2.4) the following necessary conditions : 



^2 ^ 2/xt;g - v/2(Ai + X,)vjvt - 2V2{X2 + Xs)v^ 
^ 2V2vt 

< = ^ - V-2,v, + ^-^^v^ (2.9) 

Even though, as we noted above, CP symmetry is realized at the level of the Lagrangian, 
there remains in principle the possibility for a spontaneous breakdown of this symmetry. 



^Note that the electric charge assignments for the upper and lower component fields are only conven- 
tional and can be interchanged by taking Ya = ^2, Yh = — f , entailing an exchange of the upper and lower 
components of the fermion weak doublets, without affecting the physical content. This seemingly trivial 
statement is important to keep in mind when discussing possible electric charge breaking minima of the 
potential. 
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an issue which we do not address in this paper. We can thus choose in the sequel Vd and Vt 
to be real valued; that is, we consider only CP conserving vacua for which complex valued 
Vd and/or vt can always be rotated simultaneously to real values through some unphysical 
phase redefinition of the fields. 

These equations, to which we will refer as the EWSB conditions, ensure that the vacuum 
corresponds to an extremum of the potential, [that is dV/drii\A={A),H={H) = for each of the 
ten real-valued field components denoted here by rji, {i = 1, ...10)], but one would still need 
to check that this extremum is indeed a stable, albeit local, minimum. The corresponding 
extra conditions are nothing else but the absence of tachyonic modes in the Higgs sector, 
to be considered in a later section. We just anticipate here that the latter conditions will 
enforce the signs of /i and Vt to be identical. We can thus choose in the sequel Vt > 0, 
/i > without loss of generality. Furthermore, the two free parameters m'jj and can 



now be traded for Vd and Vt through Eqs. (2.8, 2.9). In the rest of the paper we will take 



the eight parameters of the potential as being A, Xi,{i = l,.--4), fi^Vd and vf, requiring 
the correct electroweak scale will put the further constraint v = Ivl = 246GeV on 

Vd,Vt, reducing this set of free parameters down to seven. 

Let us also note that the above EWSB conditions will not necessarily imply that the 
gauge symmetric vacuum (i.e. at rji = 0) is unstable. Indeed the latter instability requires 



that < and/or mfj > which are not guaranteed by Eqs. (2.8, 2.9). Even more 
so, regimes with large fi will lead through the EWSB conditions to a very narrow gauge 
symmetric local minimum(!) so that metastability issues might have to be considered. 
[More comments about the structure of the vacua of the model will be deferred to section 

On the other edge of the spectrum, very small values of /i could be favored if one requires 



the lepton number not to be strongly violated. Indeed, the fx term in Eq.(2.4) is the only 
source of lepton number violation at the Lagrangian level and before spontaneous EWSB. 
If this term is absent the Yukawa term Eq. ( |2.5[ ) together with the other standard Yukawa 
terms imply a conserved lepton number (where the A and H Higgs fields carry respectively 
the lepton numbers I a = —2 and Zj^ = 0)j^ Then, from the lepton number assignment for 
H and A it follows that the n term violates lepton number by two units. However, this 
violation is soft since the /^-induced lepton number violating processes (corresponding either 
to loop suppressed 2—7-2 processes or to propagator suppressed multi-particle processes) 
will have to involve both the standard and neutrino Yukawa couplings . These features 
suggest that if the two seemingly independent sources of lepton number violation, namely 
the A VEV and /i, are assumed to have a common origin such as some spontaneous 
symmetry breaking of an underlying fiavor theory, then it is natural to expect /i = 0{vt) 
up to possible Yukawa coupling factors. 

2.2 Higgs masses and mixing angles 

The 10 X 10 squared mass matrix 

M'=l^\A=(A),H=iH) (2.10) 



can be recast, using Eqs. (2. 8, 2.9), in a block diagonal form of one doubly-degenerate 



The processes mediated by Eq. (2.5) and involving Higgs triplet decay or exchange are sometimes 
misleadingly dubbed 'lepton number violating'. One can check that the net overall lepton number of any 
process, comprising such decays or exchange, is conserved. This global symmetry will be violated only 



spontaneously when A acquires a VEV, that is when the Majorana mass is induced from (2.5 1 
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eigenvalue m^±± and four 2x2 matrices denoted in the following by Ai"]., ^Al■p^^^^ and 

2 

CP odd ' 



Mass of the doubly-charged field: 

The double eigenvalue m^±± corresponds to the doubly charged eigenstate and could 
also be obtained directly by collecting all the coefficients of 5^^5 in the potential. It 
reads 



m^±± (2.11) 

From now on we will denote the doubly charged mass eigenstates by if^^. 

Mass of the singly-charged field: 

The mass-squared matrix for the singly charged field is: 



This matrix is diagonalized by the following matrix TZ^' , given by : 

\ sm p cos p / 

where (3' is a rotation angle. Among the two eigenvalues of Ai"^, one is zero and corresponds 
to the charged Goldstone boson while the other corresponds to the singly charged Higgs 
boson and is given by 

, _ (t;g + 2.;,^)[2v/2/i-A4t;,] 
m^± - [Z.L6) 

The mass eigenstates and are rotated from the Lagrangian fields 0^, 6^ and defined 

by 

G^ = cos /3'(f)^ + sin /3' 6^ (2.14) 
= -sin f3'(f)^ + cos (3' (2.15) 

The diagonalization of A^^ leads to the following relations involving the rotation angle 



(2.16) 
(2.17) 
(2.18) 



7)2 


2 ^ 


= cos^/3'M^± 




2 ^ 


= ^'"f Mi. 




2 ^ 


= sin2/3'M^± 



These equat ions l ead to a unique solution for sin , cos /?' up to a global sign ambiguity. 
Indeed, Eq. ( 2.16[ ) implies V^fi— > in order not have a tachyonic state and given 
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our convention of Vt > 0. Then it follows from Eq. (2.17) that sin/5' and cos/5' should have 
the same sign. One finds 



sm/3 = , cos/3 = (2.19) 



with a sign freedom e^' = ±1, and 



tan/3' = v^— (2.20) 

Vd 



Mass of the neutral fields: 

The neutral scalar and pseudoscalar mass matrices read: 



where 



(2.22) 



2 2t;, 

These symmetric matrices are diagonalized by the following two orthogonal matrices : 

7^„=^°^"-^^^"') and n,= h'i''\^\ (2.23) 
\^ sm a cos a J \^ sm p cos p J 

where a, (3 denote the rotation angles respectively in the CVeven and CVodd sectors]^ Upon 
diagonalization of -McPe„en obtains two massive even-parity physical states /i° and 
defined by 



/^o = Cah + Saf (2.24) 
H° = -s^h + Ca^° (2.25) 

where h and are the real parts of the 0° and 6^ fields shifted by their VEV values, 

= ^{vd + h + iZi) and 5° = -^(t;* + e° + 2^2). (2.26) 
The masses are given by the eigenvalues of A^cPei,en follows, 

ml, = ^[A + C- ViA - C)2 + 4^2] (2.27) 
^?fo = ^ [A + C + v/(^ - C^)2 + 452] (2.28) 
''Hereafter, we will use the shorthand notations, Sx = s'mx and Cx = cosx, for all three angles a, /3, /3'. 
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so that niHO > rriho. Note that the hghter state h is not necessarily the hghtest of the 
Higgs sector (see section [T]). 

On the other hand, ■M.'cVodd ^^^^^ massive physical state and one massless 

Goldstone boson defined by: 



(2.29) 
(2.30) 



with masses 



m 



2 /^(^d + 4t^t ) 



A 



(2.31) 



Knowing the above eigenmasses one can then determine the rotation angles a and /3, 
which control the field content of the physical states, from the following diagonalization 
conditions: 



1. CV, 



ri 2 2,22 



B 



sin 2a; 



(m^o - m^o) 



A 2 2,22 



(2.32) 
(2.33) 
(2.34) 



2. CV 



odd- 



2 2 
sin 2/3 



f^'^d 2 2 



(2.35) 
(2.36) 

(2.37) 



Of course Eq. (2.32) should be equivalent to Eq. (2.34) upon use of + 



Eqs. (2.27, 2.28), and similarly for Eqs. (2.35) and 



, = 1 and 
S/3, C/3 will 

all be determined up to a global sign. There is however a difference between the two sectors. 



2.37). Furthermore, 



In the CVodd sector and Cfs must have the same sign as can be seen from Eq.(2.36) and 
the fact that yU > (the latter being due to the absence of tachyonic A" state,Eq.(2.35)). 
One then obtains unambiguously 



tan /3 



2vt 

Vd 



and tan 2/3 



"ivtVd 
w3 - 4f ? 



from Eqs. (2.35, 2.37), and 



2v, 



sp = ep- 



cp = ep 



Vd 



(2.38) 



(2.39) 



with a sign freedom ep = ±1. 

In contrast, the relative sign between Sa and Cq, in the CVeven sector depends on the 
values of /i as can be seen from Eqs. (2.33, 2.22). While they will have the same sign and 



tana > for most of the allowed /i and Ai, A4 ranges, there will be a small but interesting 
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domain of small /i values and tana < which we discuss in detail in section [7| One obtains 
from Eqs. ([232l-[2^ 



V2 V(A-C)2 + 452 
where eo = ±1 and e = sign[i?], and 

2B 



(2.40) 
(2.41) 



tan 2a 



A-C 



(2.42) 



Let us finally note that the angles (3 and (3 are correlated since they depend exclusively 
on Vd and Vt- For instance one always has 



tan/3 = y2tan/3' 



(2.43) 



as can be seen from Eqs. (2.20, 2.38). 



Lagrangian parameters from physical masses and couplings: 

The full experimental determination of the DTHM would require not only evidence for the 
neutral and (doubly) charged Higgs states but also the experimental determination of the 
masses and couplings of these states among themselves as well as to the gauge and matter 
sectors of the model. Crucial tests would then be driven by the predicted correlations 
among these measurable quantities. For instance one can easily express the Lagrangian 
parameters n and the A's in terms of the physical Higgs masses and the mixing angle a as 



well as the VEV's Vd,Vt, using equations (2.31, 2.13, 2.11) and (2.32 - 2.34). One finds 



Ai 



f 3 + 4:V 



1 



m 



hO 



.m 4 + 



+ cim 



2-^H± + 



sin 2a 



(m 



HO 



+ 



r.rriA - 



m 



HO) 



.171 



H± 



+ m 



H±±- 



+ 4^2""^ ^ 



.m , 



H± 



— m 



H±±. 



+ 2vt 



2 • 



r 2 2 I 2 2 1 



x/2^ 



Vt 



+ 4n 



(2.44) 
(2.45) 
(2.46) 
(2.47) 
(2.48) 

(2.49) 



The remaining two Lagrangian parameters mfj and are then related to the physical 
parameters through the EWSB conditions Eqs. (2.8, 2.9). To complete the determination 
in terms of physical quantities one should further extract the mixing angle a from the 
measurement of some couplings (see also Appendix C) and Vd and Vt from the W (or Z) 
masses. Using equations (3.53, 3.54, 2.38) one finds 
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or 



1 sin^ dw 



1 sin 2ew M| 

l + tan^/?) 2TTaQED 



vl = 



tan2/3 sin^ Ow 
(1 + I tan2 /3) AnaQED 



(2.50) 



(2.51) 



or 



sin^ 9\, 



sin^ e 



w 



21X01. 



QED 



w) 



(2.52) 



Using any of the above equations to substitute for Vd-,Vt in Eqs. (2.44, 2.49) allows us to 



obtain the Lagrangian parameters solely in terms of experimentally measurable quantities. 



Although Eqs. (2.50)-(2.52) are theoretically trivially equivalent, they involve different sets 



of experimental observables and can thus lead to non equivalent reconstruction strategies 
depending on the achieved accuracies i n the measurement of these observables. Similarly, 
trading tan/3 for tan/3 through Eq. (2.43) can be useful, depending on which of the 



two quantities is experimentally better determined through some coupling measurements. 



We should also note that Eqs. ( 2.44 )-( 2.52 ) not only allow to reconstruct the Lagrangian 
parameters from the measurable Higgs masses, a, /3, Mz and/or M^y, but can also serve 
as consistency checks among observable quantities for the model when the A's and are 
determined independently through the measurement of couplings in the purely Higgs sector 



(see also Appendix C). Finally, as can be seen from Eq.(2.52), the magnitude of Vt entails 
the deviation of the p parameter from its SM tree-level value, a point we will discuss further 
in the following section. 



3 Miscellaneous constraints 

3.1 Constraints from electroweak precision measurements 

In the Standard Model the custodial symmetry ensures that the p parameter, p = — 



cos^ dw 

is equal to 1 at tree level. In the DTHM one obtains the Z and W gauge boson masses 



readily from Eq.(2.7) and the kinetic terms in Eq.(2.1) as 



^ (j!±m±M) = 9!M±i!^) (3 53) 

^ 4 4cos2^w' ^ ■ ' 

Ml = ^!(!i±M) (3.54) 



whence the modified form of the p parameter: 



vl + H 

vl + ^vl 



P = ^ 1 (3-55) 



and actually p < 1 at the tree-level. Since we are interested in the limit Vt ^ Vd we rewrite 
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p~l-24 = l + '5p (3.56) 

with 5p = —2\ < and a/w J + 2v^ = 246 GeV. Thus the model will remain viable 
as far as the experimentally driven values of 6p are compatible with a negative number. 
The implication for the DTHM has already been studied in the literature [17]. Here we 
only discuss briefly this point taking into account the latest updates of the electroweak 
observables fits as reported by the PDG |tl9j. One should compare the theoretical value 
with the experimental value after having subtracted from the latter all the known standard 
model contributions to the p parameter. The quoted number after this subtraction, po = 
1.0008!l^Q QQQ7 obtained from a global fit including the direct search limits on the standard 
Higgs boson, is not compatible with a negative Sp and would exclude the DTHM. However 
at the 2 a level, one obtains po = 1.0004;!^q;qq^^ [19j, which is again compatible with Sp < 0. 
Moreover, relaxing the Higgs direct limit leads to po = 1.0008^o;Qojg, again compatible with 
6p < 0. From the last two numbers one gets an upper bound on Vt of order 2.5 — 4.6 GeV. In 
the present study we will thus be contented by the conservative assumption that an upper 
bound of 2.5 GeV guarantees consistency with the experimental constraints. We should 
note, though, that the tree-level DTHM value of 6p being of order 10~^, it is legitimate 
to ask about the effects of radiative corrections to this quantity within the DTHM. As far 
as we know, radiative corrections to 6p are not available in the literature in the case of 
Ya = 2 that concerns us here, while several studies have been dedicated to this question in 
the framework of a Ya = triplet Higgs [20l EH [22l [23] , In |20] , it has been shown that the 
tree level bound on the triplet vev could be pushed to higher values by one-loop radiative 
corrections. Whether this happens also in our case is still to be investigated and deserves 
a study that is out of the scope of the present paper, including for that matter all other 
LEP/SLC SM observables. 



3.2 Absence of tachyonic modes 



From Eq. (2.31), the requirement that mj^ should be positive implies pvt > 0. The same 



positivity requirement in the singly charged and doubly charged sectors, Eqs. (2.13, 2.11), 



together with our phase convention > discussed in section [2| lead to the following 
bounds on p: 



p> 



p > 
p > 



V2 



+ V2 



(3.57) 
(3.58) 

(3.59) 



The tachyonless condition in the CVeven sector, Eqs. (2.27, 2.28 ), is somewhat more involved 
and reads 

^/2pvl + Xvlvt + 4(A2 + XzWt > (3.60) 
-Sp'^vt + V2p{Xvj + 8(Ai + X4)v^) + 4(A(A2 + A3) - (Ai + X4f)vf > (3.61) 

The first of these two equations is actually always satisfied as a consequence of Eq. (|3.57 ) 
and the boundedness from below conditions for the potential (see section|4]and Eq. (4.21 )). 
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The second equation, quadratic in /x, will lead to new constraints on in the form of an 
allowed range 



/i_ < /i < /i^ 



(3.62) 



The full expressions of fi± and a discussion of their real-valuedness are given in appendix A. 
Here we discuss their behavior in the regime Vt <^ Vd- In this case one finds a vanishingly 
small /i_ given by 



and a large given by 



((Ai + A4)'-A(A2 + A3)) 



+ 0{vl 



(3.63) 



+ V2(Ai + Xi)vt + 0{vt 



(3.64) 



Depending on the signs and magnitudes of the A's, one of the lower bounds (3.57 - 13.59[ ) 



or will overwhelm the others. Moreover, these no-tachyon bounds will have eventually 
to be amended by taking into account the existing experimental exclusion limits. This is 
straightforward for y4°, and H^^. We thus define for later reference 



max 



[m 



2^2 



+ 



m 



.V2^^ vl 



V2vt 



(m^±±) 



exp 



(3.65) 



where {mA)exp, (^H±)cxp, (^H±±)cxp denote some experimental lower exclusion limits for 
the Higgs masses. Eqs.( |3.57 - 3.59) are then replaced by 



II > /imin (3.66) 

in order for the masses to satisfy these exclusion limits. Similar modifications on fi± taking 
into account experimental exclusion limits in the CVeven sector are more involved and will 
be differed to section [7] after having established the theoretical upper (lower) bounds on 
the /i° {H^) masses. Furthermore, the upper bound /i+ will be instrumental in determining 
the maximally allowed values of the six Higgs masses ■mHO,mA,'mH±,mH±± as we will see 
in section [71 



3.3 The vacuum structure 

Obviously, violation of any of the constraints discussed in the previous subsection is a 
signal that the would-be electroweak vacuum is not a minimum (but rather a saddle-point 
or a local maximum) for the given set of values A, Xi,Vd,Vt when /i is either very small or 



very large. However, since Eqs.(2.8, 2.9) are non-linear in Vd,Vt, it could still be possible 
to find a different set of values for the same input values of m^,M^, where the 

true electroweak minimum is obtained at a lower point of the potential than the previous 
one. More generally, and depending on the values of the parameters of the potential, 
one expects on top of the electroweak minimum a rich structure of extrema that can 
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affect tlie interpretation and viability of tliis minimum and tlius possibly lead to additional 
constraints on these parameters. A complete study of such extrema can be very involved 
since the potential depends on ten independent real fields. Here we only provide a partial 
qualitative discussion. 



Upon use of Eqs.(2.7, 2.8, 2.9) in Eq.(2.4) one readily finds that the value of the 



potential at the electroweak minimum, (\^)ewsb, is given by: 

(I^)ewsb = + 4(A2 + X3)vt + Av^vtiiXi + X,)vt - v^/i)) (3.67) 

16 

Since the potential vanishes at the gauge invariant origin of the field space, Vh=o,a=o = 
0, then spontaneous electroweak symmetry breaking would be energetically disfavored if 
(V)ewsb > Oj^ One can thus require as a first approximation the naive bound on n 

/X < /i^a, ^ -A^^ + (Ai + A4)^ + 0{v^) (3.68) 
4v2 Vt V2 



to ensure I^wsb < 0. As can be seen from Eq.(3.64) one has either /imax < yU+ or /^max > /^+ 



depending on the sign of Ai + A4. But for all practical purposes /^max — yU+ in the regime 
Vt/vd ^ 1, so that the proviso stated above concerning the relevance of the tachyonless 
conditions is weakened for the upper bound which can be replaced by /imax- There is. 



however, yet another critical value of fi. As mentioned at the end of section 2.1, and 
—m'jj can be both positive for sufficiently large values of /i thus making the gauge invariant 
point if = 0, A = a local minimum. This happens when /i > where 

'"''^47l'''^^^'^^'^27f ^^'^^^ 
If A1 + A4 > then fiH < /imax < /i+- To delineate some consistency constraints in this case, 
it would be necessary to look more closely at the decay rate from a metastable gauge invari- 
ant vacuum to the EWSB vacuum, if < /t < /imax, and vice-versa, from a metastable 
EWSB vacuum to the gauge invariant vacuum when /tmax < /i < /i+. Fortunately, however, 
these configurations altogether are already excluded if we take into account the experimen- 
tal mass limits on the Standard Model Higgs. Indeed, as will be shown in sections [7] and 
[s} the lightest CVeven Higgs state /i" becomes purely standard model-like for such large 
values of /i, irrespective of the values of the couplings A, A,, while m^o becomes very small 
for these values (e.g. = a/3(Ai + X4)vt for /i = /i//) and thus experimentally excluded. 



Nonetheless, the structure of the potential Eq.(2.4) is sufficiently rich to provide danger- 
ous extrema configurations which are not excluded by the above mentioned experimental 
limits. We exhibit here, without much details, one example among a manifold of possibili- 
ties. There is an extremum in the field space direction defined by Re 0° = Re 0"'" = ^ and 

Re5° = — Re 5"'""'" = and all other fields put to zero. 
This requires 

/i = (3.70) 

ml = ^-{Xvf + {2X^ - X,X') (3.71) 
Ml = -X,vf - (2A2 + Xs)vf (3.72) 



^We should, however, keep in mind the possibihty that a long-hved metastable vacuum could still be 
physically acceptable, even when (T^)ewsb > 0, thus altering our constraints; these issues are not addressed 
further in the present paper. 
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Note that this direction, and thus the corresponding extremum, breaks spontaneously 
charge conservation. We will refer to this extremum as charge breaking (CB). Furthermore, 
in contrast with the EWSB point, Eqs.(2.8, 2.9), here fi is not a free parameter. We 



can then seek for a region in parameter space where this CB extremum coexists with an 
EWSB minimum, and check what happens at the gauge invariant extremum point as well. 



Requiring Eqs.(2.8, 2.9, 3.70 - 3.72) to be simultaneously satisfied leads to correlations 
among Vd,Vt,v^,v^. These lead in turn to constraints on the A, A, parameter space in 
order for all these vev's to be real valued (modulo gauge transformations), together with 
the immediate constraint X4V^ < originating from Eq.(3.70) and > 0|^ The ensuing 



correlations allow to write mf^o, (^)ewsb and (I^)cb (the value of the potential at the CB 
extremum) in the following form: 



{V) 



mlo 



EWSB 



(A(2A2 + A3) + 2A^ + Ai(A4 
-(A(2A2 + A3) + 2A^ + Ai(A4 - 2Ai)) x 



(A(2A2 + A3) + Ai(-2Ai + A4)) 



4AA1 



{V) 



{AX\-2X{2X, + X,)~Xl)'^ + 0{vi) 



c4 



CB 



(3.73) 



(3.74) 



(3.75) 



Various interesting conclusions can be drawn from the above equations. As can be seen from 



Eq.(3.73), a physical 



, i.e. m^o > 0, implies a positive m^j and thus an unstable gauge 
invariant point at the origin of the fields {H = 0, A = 0). Furthermore, in the consistent 
(A, Aj) domain (given in footnote 2) m^o is indeed positive and, furthermore, one finds 



from Eq.(3.74) that (I^)ewsb < 0. The EWSB vacuum is thus energetically favored over 
the gauge symmetry preserving one which lies at = 0. It then remains to compare the 



EWSB point with the CB point. Close inspection of Eq.(3.75) shows that {V)cb > in 



all the (A, Aj) domain given in footnote 2, if and only if A4 < 0, in which case the EWSB 
is energetically favored over the CB. However, if A4 > (and thus < 0), there are 
regions in the (A, Aj) consistent domain where (V^)cb < 0, provided that 4A4 < A(2A2 + A3). 
Moreover, the potential at this CB point becomes much deeper than at the EWSB point 
since we are in the regime Vt <^ This is a dangerous configuration since it makes 

the EWSB vacuum potentially very short-lived due to tunneling effects [Ml EH] . We stress 
here that this EWSB point is a true local minimum in this configuration, i.e. there are 
no tachyonic Higgs states which could have signaled its non-relevance beforehand. [This 
is easily seen from the fact that is non-tachyonic and is the lightest Higgs state when 
/i ~ \v^ \ ^ Vt] see also section [7j] Even more so, the potential is bounded from below, as 
can be shown by comparing the corresponding (A, Ai) domain with the boundedness from 
below constraints that we will derive in the following section. We have thus exhibited an 
example of a configuration where /i can be very large, consistent with the experimental 

mass limit, and a fortiori with all the non-tachyon constraints, corresponding locally to 
an acceptable EWSB vacuum, but still non-viable due to the existence of lower (charge 
breaking) points akin to what happens in two-Higgs-doublet models (see for instance [2S])- 

We end this section with a general comment concerning the neutrino mass see-saw 
mechanism. The common lore is to assume a GUT origin for fi and Ma, and taking 



^Working in the regime Vt ^ Vd, 



and keeping only terms 0{vt), the constraint in the (A, A,) 



space satisfying all these requirements is found to be Ai < ^(A4— •\/8A(2A2 
and 2X2 + A3 > for a bounded from below potential (see section 111). 



A, 



17A|). Note that A > 
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fi ~ Ma ~ C(Mgut) leads through Eq.(2.8) naturally to a tiny vt- However, as noted 
in the introduction we do not commit in the present study to specific high energy physics 
scenarios, so that Ma and/or fi could be smaller than a hypothetical GUT scale. It is 
then interesting to note that even in this case a kind of see-saw is actually still at work 
model-independently due to the dynamics of the potential. This is simply due to the form 
of the fi upper bound /i^, Eq.(3.64): the larger is /i the smaller should be Vt in order to 
avoid a tachyonic h^. For instance, taking A ~ 0.5 and /i+ ^ 2 x lO^^GeV leads to Vt — leV 
and Ma ~ lO^^GeV. 



4 Boundedness of the potential 

A necessary condition for the stability of the vacuum comes from requiring that the poten- 



tial given in Eq. (2.4) be bounded from below when the scalar fields become large in any 
direction of the field space. The constraints ensuring boundedness from below (BFB) of 
the DTHM potential have been studied in the literature so far only partially (see e.g. |18j). 
and at the tree-level. It would thus be somewhat premature to invoke possible quantum 
modifications of these constraints before fully settling first the tree-level issue. This section 
is devoted to this issue and aims at deriving, at the tree-level, the complete necessary and 
sufEcient BFB conditions valid for all directions in field spacej^ 



Obviously, at large field values the potential Eq. (2.4) is generically dominated by the 



part containing the terms that are quartic in the fields. 



V^^\H, A) = ^{H^Hf + Xi{H^H)Tr{A^A) + X2{TrA^Af + X^TriA^AY + X^H^AA^H 

(4.1) 

The study of V^'^\H, A) will thus be sufficient to obtain the main constraints. To obtain 
BFB conditions it is common in the literature to pick up specific field directions or to put 
some of the couplings to zero. Consider for instance the two following cases: 

1. in the absence of any coupling between doublet and triplet Higgs bosons, i.e. Ai = 
A4 = 0, it is obvious that 



A>0 & A2>0 & A3>0 (4.2) 

will ensure that the potential is bounded from below. 

2. if one picks up the field space directions where only the electrically neutral compo- 
nents are non vanishing, one finds 

Vo^'^ = ^|0r + (A2 + A3)|5r + (Ai + A4)|0T|5T (4.3) 



In order for the potential to be bounded from below in this sub-space, V^*^^"* should be 
positive for any values of and \6^\ including when one or the other is vanishing. 
The latter cases imply the necessary conditions A > and A2 + A3 > 0. It is then 



possible to rewrite Eq.(4.3) in the form 



^We will thus not address in this paper the possibility that loop corrections could lift the potential in 
some otherwise unbounded from below directions, nor the issues related to metastability of the vacuum 



which could relax some of the constraints. See also section 3.3 
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(4) 



Since the first term is non-negative and vanishes in the direction 

1 00 1 2^ 1^0 1 2 _ 2a/(A2 + A3)/A, then the necessary and sufficient conditions for the 

BFB of the potential in this direction are 



A > 
A2 + A3 > 
Ai + A4 + VA(A2 + A3) > 



(4.5) 



As it will become clear later on in this section, the conditions in case 1 are sufficient but 
not necessary, even for this special case. Furthermore, while the conditions in case 2 are 
necessary and sufficient for the corresponding direction, they obviously remain necessary 
for the general potential but it is a priori not clear whether they can be sufficient. By 
looking at other special directions in 2-field and 3-field directions we will show that they 
are generally not sufficient. Before doing so, let us first point out a more convenient method 



to obtain positivity constraints like Eq.(4.5) directly from Eq.(4.3) rather than writing it 



first in the form of Eq.(4.4). The potential Eq.(4.3) can be cast in the form 



V{x) = a + bx^ + cx' 



(4.6) 



by the change of variable x = Since x is by definition real valued and the moduli 

10*^1 and I (5° I can have any value, then the problem of finding the necessary and sufficient 



BFB conditions for Eq.(4.3) is equivalent to finding the conditions on a,b,c such that 



V{x) > for any x ^ [0, oo[. Since V{x) has no linear or cubic terms in x it is easy to 
find these conditions by studying V{x) as a bi-quadratic function: 



a > 
c> 
b + 2y/ac > 



(4.7) 



Applied to Eq.(4.3) these conditions reproduce immediately Eq. (4.5). We can now easily 



study other field directions. For instance, the direction where only 6~^~^ and 
vanishing yields 



^ = (A2 + A3)|(5++|^ + Ai|5 



++|2 I ,0|2 



A 



/.0|4 



for which the BFB constraints are readily obtained from Eq.(4.7) as 



A>0 & A2 + A3>0 & Ai + v^A(A2 + A3) > 0. 
Similarly, if we consider the field direction with non-vanishing 5"*" and 



V = {X2 + ^)\S^\' + {Xi + ^)\6^\'\<P'-\' + 



A 



then 



are non 



(4. 



(4.9) 



(4.10) 



and the corresponding BFB conditions read 
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A > & Ao 



> 0. 



(4.11) 



^>0& A, + ^ + yA(A. + ^ 

It is then obvious that these two sets of conditions are neither equivalent nor contained 
in the conditions of Eq.(4.5). This shows that the BFB conditions derived only from the 



neutral direction Eq.(4.3) are neither necessary nor sufficient to ensure boundedness from 



below of the full potential Eq.(2.4). In Appendix B we have listed the potentials for all the 



field directions with only two non-vanishing fields together with the corresponding BFB 
conditions. Adding these conditions we come closer to the real sufficient and necessary 
conditions. But one can get more conditions by going now to field directions where 3 fields 
are non-vanishing. We give the exhaustive list of all these 3-field directions potentials in 
Appendix B. In these more complicated configurations, an iteration of the method described 
above allowed to treat all of them, although the results become somewhat complicated and 
not so telling. For instance the 3-field direction with non-vanishing 0°, 0+, 6~^, see Eq.(B.24) 
yields some of the simplest BFB conditions 



A > A 2A2 + As > A v/A(4A2 + 2A3) + 2Ai + A4 > A 
(2X{2\2 + A3) > (2Ai + V 2Ai + A4 > 0) 



(4.12) 



where A, V stand respectively for the logical AND, OR. These conditions are obtained by 
first defining the reduced variables Xi = \4>~^\/\(p'^\,X2 = |<^^ 1/10° 1 5 and then using iteratively 



the constraints Eqs. (4.7). By the same method we could obtain even more complicated 
BFB conditions given in Eqs. (B.26 - B.35). Analyzing them numerically we confirm 
that Eqs. (4.5) are far from being the full story. However, and despite their apparently 



complicated structure the intersection of the regions they delineate in the space of the 
A's has a form similar to equations (4.5) and (B.14). Moreover, the true BFB conditions 
will be obtained only if all field directions are taken into account, up to some arbitrary 
SU (2) X U (1) gauge transformations, but in this case the method used so far is not tractable 
anymore. 

To proceed to the most general case, we adopt a different parameterization of the fields 
that will turn out to be particularly convenient to entirely solve the problem. Without loss 
of generality we can define: 



Tr{A^A) 
Tr(AtA)V(TrAtA)' 
{H^AA^H) /{H^HTrA^A) 



^WH + Tr AtA 

2 2 
r cos 7 

2 • 2 
r sm 7 

c 



(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 



(where we adopted here a parameterization similar to the one used in [21] to study two- 
Higgs-doublet models, although for the latter models the problem is not fully solved by 
such a parameterization). Obviously, when H and A scan all the field space, the radius r 
scans the domain [0, oo[ and the angle 7 G [0, |]. Moreover, one can show that 



0<^<land-<C<l 



(4.1^ 
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With this parameterization it is straightforward to cast V^^^ (if, A) in the following simple 
form, 



l^(4)(r,tan7,e,C) 



4(l + tan2 7)2 



(A + 4(Ai + eA4) tan^ 7 + 4(A2 + CA3) tan^ 7) (4.19) 



Due to the bi-quadratic dependence in tan 7, one can indeed consider only the range 
< tan 7 < +00 in accor dan ce with the above stated range for 7. We have thus written 
V^'^'^ in the form of Eq. (4.6). Boundedness from below is then equivalent to requiring 



V^'^^ > for all tan 7 E [0, oo[ and all ^, ( satisfying Eq.(4.18). Now applying directly the 



conditions Eq.(4.7) one obtains 



A > & A2 + CA3 > & Ai + eA4+ v/A(A2 + CA3) > VC e 1], Ve[0, 1] (4.20) 

Due to the monotonic dependence in ( and ^, it is easy to show that these conditions 
can be rewritten as. 



As 

A>0 & A2 + A3>0 & A2 + — >0 



& Ai + VA(A2 + A3) > & Ai + ^/ A(A2 + y ) > 



& Ai + A4 + VA(A2 + A3) > & Ai + A4 + ^ A(A2 + y ) > (4.23) 

We stress here that the above conditions ensure BFB for all possible directions in 
field space and thus provide the most general all directions necessary and sufEcient BFB 
conditions that solve completely the issue at the tree-level. Note that all the 2-field direc- 



(4.21) 
(4.22) 



tions conditions given in Eqs.(B.ll - B.15 ) are special cases of the above conditions. We 
also checked numerically that this is the case for all the ten 3-field directions conditions 



Eqs. (B.26 - B.35). 



5 Unitarity constraints 

Constraints on the scalar potential parameters can be obtained by demanding that tree- 
level unitarity be preserved in a variety of scattering processes: scalar-scalar scattering, 
gauge boson-gauge boson scattering and scalar-gauge boson scattering as was initially done 
for the SM [281 129| 130] . The generalizations of such constraints to various extended Higgs 
sector scenarios have been studied in the hterature, see for instance [311 [331 IMj- Here 
we treat such constraints in the DTHM at the tree-level, limiting ourselves to 2-body scalar 
scattering processes dominated by quartic interactions. This is justified by the fact that 
we are interested in the leading unitarity constraints, that is in the limit where ^/s is much 
larger than any other mass scale involved. In particular this means that we disregard 
here unitarity constraints that would involve the /i parameter when the latter is very 
large. Indeed, this parameter contributes to the scalar scattering processes through the 
cubic interactions entering the Feynman diagrams with scalar exchange in the s, t and u 
channels. Furthermore, the ratio fi/vt controls the size of the exchanged scalar masses so 
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that some of the aforementioned diagrams can be important in the vicinity of the resonance 
pole in the hmit of large a/s ~ fiVd/vt. 

In order to derive the unitarity constraints on the scalar masses we adopt the basis 
of unrotated states, corresponding to the fields before electroweak symmetry breaking. 
The quartic scalar vertices have in this case a much simpler form than the complicated 
functions of Aj, a and (3 obtained in the physical basis {H^^, H^, G^, h^, H^, and 
G^) of mass eigenstate fields. The S'-matrix for the physical fields is related by a unitary 
transformation to the S'-matrix for the unrotated fields. Close inspection shows that the 
full set of 2-body scalar scattering processes leads to a 35 x 35 S'-matrix which can be 
decomposed into 7 block submatrices corresponding to mutually unmixed sets of channels 
with definite charge and CP states. One has the following submatrix dimensions, structured 
in terms of net electric charge in the initial/final states: S^^\6 x 6), S^'^\7 x 7) and 
S^^\2 X 2), corresponding to 0-charge channels, S'*^^^(10 x 10) corresponding to the 1-charge 
channels, S'*^^''(7 x 7) corresponding to the 2-charge channels, S'''^''(2 x 2) corresponding to 
the 3-charge channels and S'*^''^(l x 1) corresponding to the unique 4-charge channel. The 
corresponding T-matrix submatrices T^^\ ...,T^'^^ -with a momentum conservation factor 
(27r)^5^(^ momenta) properly factored out- are then easily extracted using the pure scalar 
quartic interactions expressed in terms of the non-physical fields 0^, 6^, S^^, h, C,^, and 
Zj(i=l,2) as hsted in Appendix C. 

One can then in principle extract the unitarity constraints on each component of the 
T-matrix through the unitarity equation which we write here in a shorthand form as 



(5.1) 



where J denotes symbolically the phase space integral over each intermediate state channel 
(see for instance [35]). However, it proves more efficient to define a modified matrix in 



such a way that its diagonalized form still satisfies Eq. (5.1). The usual unitarity bound 
on partial-wave amplitudes that is valid for elastic scattering would then apply readily 
to all the eigenvalues, thus encoding indirectly the bounds on all the components of the 
T-matrix|^ The proper redefinition is a T-matrix having the same entries as T but with 
an extra 1 / a/2 factor for each initial or final state channel having two identical particles. 



T now satisfies Eq. ( |5.1 ) with the same phase space for all channels and a true matrix 



multiplication of T by T^. Its diagonalized form will thus satisfy the same equation. 
Defining = iT^^\ with n = 1, ■ ■ ■ , 7, we give hereafter the resulting submatrices whose 
entries correspond to the quartic couplings that mediate the 2 — )■ 2 scalar processes. These 
submatrices are hermitian, thus the sought for eigenvalues will all be real-valued. 

The first submatrix TWi corresponds to scattering whose initial and final states are one 
of the following: (0+5", 5+0^, Z1Z2, h^^). With the help of Appendix C one 

finds, 



^This, however, cannot be achieved in general by simply diagonalizing T, since on the right-hand side 



of Eq. (5.1) the phase space factor is not the same for all the 2-particle channels, even in the high energy 



massless limit we are considering. It picks up a factor 1/2 only for internal states with identical particles 



so as to avoid double counting. The right-hand side of Eq. ( 5.1 1 is thus not a proper matrix multiplication 



of T by , a fact emphasized by the quotation marks in the equation 
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(5.2) 



where A,^ = Aj ± Xj. We find that M.i has the following 3 double eigenvalues: 



ei 

62 



Ai + A4 

Ai 

1 



63 = -(2A1 + 3A4) 



(5.3) 
(5.4) 

(5.5) 



The second submatrix M.2 corresponds to scattering with one of the following initial 
and final states: (0+0-, 5+5", 5++ 5—, ^, S^), where the ^2 accounts for 

identical particle statistics. Again, with the help of Appendix C, one finds that M.2 is 
given by: 



M2 = 
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(5.6) 



where A14 = 2Ai + A4 and A23 = 2A2 + A3. Despite its apparently complicated structure, 
one can easily determine the seven eigenvalues of AI2 as follows: 
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A 



/2 — 2A2 

/3 = 2(A2 + A3) 

1 

4 



a± 



[A + 4A2 + 8A3 ± \/ (A - 4A2 - 8A3) 



16A2 



M3^ 



^ 

2A+3, 



(5.7) 

(5.8) 
(5.9) 

(5.10) 



= |[3A + I6A2 + I2A3 ± V(3A - I6A2 - 12A3)2 + 24(2Ai + X^^ ] (5.11) 
The third submatrix M.^ corresponds to the basis {hZ\^^^ Z2) and is given by: 

A 



(5.12) 



with eigenvalues k\ — j\ and k^ — fz- 

The 1-charge channels occur for two-by-two body scattering between the 10 charged 
states (/i0+, CV+, ^10+, ^20+, ^K^^, ^2^+, 5++0-). The 10x10 

submatrix AI4 obtained from the above scattering processes is given by: 
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iA4 


2v^ 


A4 


2 


— iA4 


2^2 


— iA4 


2^ 


A4 










2 





2A14 


2v^ 






2A^3 








-y2A3 








^A 

















2A^3 


-^V2A3 





-^A3 





iv^A3 


2AJ3 


-A4 

2 





iA4 

~Y 














-A4 

2 


-iA4 

~Y~ 




Ar4 / 



(5.13) 



As one can see, this matrix contains many vanishing elements, and the 10 eigenvalues are 
straightforward to obtain analytically. They read as follows: 
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di = ei 


(5.14) 


d2 — 62 (twice) 


(5.15) 


ds = 63 


(5.16) 


d4 = fi 


(5.17) 


dr — fo 


I tJ. -LO I 


"6 — J3 


yo.i-ij J 




(5.20) 




(5.21) 



The fifth submatrix Ai^ corresponds to scattering with initial and final states being one 
of the following 7 sates: (^,^,5+0+,^++^°, 6++Z^, 6++h). It reads, 







2 

A23 




-A3 

iA3 










Al4 




iX^ 

\ 




-A3 -iA3 





2A2 









2A2 









-iXi 

~Y~ 




Ai 







-A4 

2 





Ai 



(5.22) 



Cl 


= ei 


C2 


= 62 


C3 


= /l 


C4 


= /2 


C5 


= /3 


C6 


= ^7 


C7 


= 2A2 — A3 



and possesses the following 7 distinct eigenvalues: 

(5.23) 
(5.24) 
(5.25) 
(5.26) 
(5.27) 
(5.28) 
(5.29) 

There are also triply-charged states.The submatrix M.^ corresponding to this case generates 
the scattering with initial and final states being one of the following (5++0+,5++5+), and 
is given by : 

with eigenvalues k\ = ci and /c2 = fs- Finally, it is easy to check that there is just one 
quadruply-charged state -^^^'^^'^'^ leading to 

= h (5.31) 
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with /s as eigenvalue. 

From the usual expansion in terms of partial- wave amplitudes aj, we write, following 
our notations, 

A^C^/) =ifkf = 16m5]](2J+ l)ay^Hs)Pj(cos^) (5.32) 

J>0 



where Ai^''-^^ denotes the entries of the Ai matrix, the subscripts k and / run over all possi- 
ble initial and final states of the above 35-state basis, 6 denotes the scattering angle of the 
corresponding processes and the Pj's are the Legendre polynomials. Since we considered 
only the leading high energy (massless limit) contributions that are s and 6 independent, 
all the partial waves with J ^ vanish and one is left with 



16n 



for each channel. The S-matrix unitarity constraint for elastic scattering |a, 



ter natively |-Re(aQ 



{kk)i 
I 



(5.33) 



< 1 (or al- 



< 



[37] ) applies to the diagonal entries of . It encodes as well 



the constraints for non-elastic scattering, provided that it is applied to the eigenchannels of 



the 35-state basis as noted previously. Thus, this constraint translates through Eq. (5.33) 



directly to all the eigenvalues we determined above. We defer to the next section, Eqs. (6.4 



6.13) the list of all the resulting constraints. 



6 Combined unitarity and potential stability constraints 

Let us first recall all the constraints obtained in sections H] and |5l 



BFB: 



A>0 & A2 + A3>0 & A2 + -r>0 



& Ai + VA(A2 + A3) > & Ai + ^/a(A2 + y) > 



& Ai + A4 + VMMTUj > & Ai + A4 + \/ A(A2 + y ) > 



unitarity: 



Al + A4I < KTC 
All < KIT 

2Ai + 3A4I < 2Kn 
A| < 2k7v 

K, 

21-2 

K 

A2 + A3 < -71 



A, 



A + 4A2 + 8A3 ± -^(A - 4A2 - 8A3)2 + 16A| I < Aktt 

3A + I6A2 + I2A3 ± V(3A - I6A2 - 12A3)2 + 24(2Ai + A4)2 | < Ann 
2Ai - A4I < 2k7t 

2X2 — A3I < KTT 



(6.1) 
(6.2) 
(6.3) 



(6.4) 
(6.5) 
(6.6) 
(6.7) 

(6.8) 
(6.9) 

(6.10) 

(6.11) 
(6.12) 
(6.13) 
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where we introduced the parameter k that takes the values k = 16 or 8, depending on 
whether we choose |ao| < 1 or |_Re(ao)| < | as pointed out at the end of section [sj 
Working out analytically these two sets of BFB and unitarity constraints, one can reduce 
them to a more compact system where the allowed ranges for the A's are easily identified. 
One can obtain a necessary domain for A,A2,A3 that does not depend on Ai and A4, by 
considering simultaneously Eqs.(6.7 - 6.13) together with Eq. (6.1). It then turns out that 



Eqs.(6.8, 6.9) as well as the lower part of Eq.(6.13) are weaker than the actually allowed 



domains for A2, A3, and similarly Eq.(6.7) is weaker than the constraint on A coming from 



Eq.(6.11). We find. 



< A < -KTT 

- - 3 



A3 



A2 + A3 > & A2 + ^ > 



A2 + 2A3 < -vr 



K 



4A2 + 3A3 < -vr 
2A9 — A3 < /tvr 



(6.14) 

(6.15) 

(6.16) 

(6.17) 
(6.18) 



We stress here that the above constraints define the largest possible domain for A, A2, A3 for 
any set of allowed values of Ai, A4, although Eqs.(6.10, 6.11) have been used to determine 



this domain. It is noteworthy that the upper bound on A Eq.(6.14) is reduced by a factor 
3 with respect to the naive expectation Eq. (6.7). Studying further Eqs.(6.10, 6.11), one 



can rewrite them in the following simple form where the dependence on Ai,A4 has been 
explicitly separated from that on A, A2, A3: 



K 



IA4I < min J (A ± 2K7r)(A2 + 2A3 ± -vr) 



|2Ai + A4I < \l2{\ - ^/€vr)(4A2 + 3A3 - ^vr) 



(6.19) 
(6.20) 



where Eqs.(6.14, 6.18) have been used in deriving Eq.(6.20)^ Various comments are in 
order here. First, to obtain the full domain for Ai, A4, one has to add to Eqs.(6.19, 6.20) 



equations (6.4 - 6.6, 6.12) as well as Eqs.(6.2, 6.3). Thus for each set of values of A, A2, A 



the allowed domain for Ai, A4 is easily determined as the overlap of a set of linear bands, 
as illustrated in Fig. 1. 

As stated earlier, Eqs. (6.15 - 6.18) define the largest possible domain for A2, A3 allowed 



by the combined unitarity and BFB constraints. The reason is seen from Eqs. (6.19, 6.20) 
that are the only extra constraints on A2,A3 depending on the actual values of A, Ai,A4. 
As one can easily check, these constraints become trivially satisfied when Ai = A4 = and 
thus correspond to the case of largest A2, A3 domain. For each set of non vanishing values 



for Ai, A4 the domain of A2, A3 given by Eqs. ( |6.15 - 6.18) will be further reduced according 
to Eqs. (6.19, 6.20). We illustrate the largest (A2, A3) domain in Fig. 2. 



In writing Eq. (6.20) we relied on the fact that the minimum of i/2(A± 4K7r)(4A2 + 3A3 ± |7r) 



is given by \J'^{^ — |K7r)(4A2 + 3A3 — ^tt) in all the domain allowed by A, A2,A3. In contrast, 



mm ^(A±2«;7r)(A2 + 2A3± f tt) appearing in Eq. te.l9b cannot be written in a closed form in this domain 
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Figure 1: An illustration of a section of the (Ai,A4) domain (blue) in units of k as 
determined by Eqs. (6.4 - 6.6, 6.12, 6.19, 6.20), where we fixed A = A2 = A3 = 0. Adding 



the BFB constraints Eqs. (6.2, 6.3) one obtains the reduced domain shown (light purple). 



To summarize, the boundaries of the combined unitarity and general BFB domains for 



the five couplings are now given by the reduced set of equations (6.2 - 6.6, 6.14- 6.20) 



which moreover have an analytically simpler form. In particular, one readily finds from 
Eq. (6.20) that saturating the unitarity bound on A, i.e. A = |K7r, reduces the two- 
dimensional (Ai,A4) domain to the one-dimensional (straight line) A4 = — 2Ai. This, as 
well as other features, will be useful in determining lower and upper bounds on the Higgs 
masses in the next section. 



7 Higgs mass theoretical bounds 

In this section we rely on the results of the previous sections to study the theoretically 
allowed ranges of the Higgs masses when varying the A'^s and the n parameter in their 
allowed domains. Rather than assuming that fi is very large, i.e. of the order of the GUT 
scale together with fi ^ Ma, we will study all the phenomenologically allowed range. We 
stress here that even very small values of /i are consistent with a tiny value of Vt necessary 
for realistic neutrino masses (and 0{1) Yukawa couplings) provided that we take into 
account consistently Eq. ( 2.8[ ). 

Let us first describe qualitatively the generic behavior of the masses when fi is varied. 
We will show that, as a function of /i, the mass features a maximum m™j^^ for a specific 
value II = fic- This maximum will translate into an upper bound on m/jO when the unitarity 
bound on A is saturated. Similarly, the mass reaches a minimum m^™ at a nearby value 
which we momentarily denote also fi = fic for the sake of the qualitative discussion. In 



25 



1.0 



0.5 - 



St 

o 

0.0 

c 

3 

C 



-0.5 



-1.0 




-0.5 



0.0 0.5 

A2 in units of k 



1.0 



Figure 2: We illustrate here the largest A2 — A3 domain allowed by the combined unitarity 
and BFB constraints in units of n. This domain corresponds to Eqs. (6.15 - 6.18) and is 
attained for Ai = A4 = in which case Eqs. (6.19, 6.20) are trivially satisfied. As can be 



seen from Eqs. (6.19, 6.20), a smaller domain obtains as soon as Ai and/or A4 are non-zero, 



irrespecive of the value of A satisfying Eq. (6.14) 



the range /i < /Zc, is the heaviest among all the Higgses, decreasing very slowly with 
increasing /x towards its minimum value m™o", while rriho increases very quickly with /i to 
The other Higgs masses can have various hierarchies and in particular the unusual 



one where the mH±± is the lightest state, mH±± < rnH± < — rnho. In contrast, in the 
range fic < < /imax, now increases quickly with while rrihO decreases very slowly 
from its maximal value. This sharply different behavior of rrifiO and ttiho below and above 
/ic can be traced back to the smallness of Vt- We illustrate numerically such a behavior in 
Fig. [3] where the seemingly constant m^o for fi > jj^c and constant m^o for jj, < jj^c is an 



artifact of the very small ratio Vt/vd- In fact. 



m 



hi. 



cf. section 



3.2 



of /ic and reaches zero when /i 

off at /i = /imin decreasing very slowly until its minimum value at /i 



is decreasing very slowly to the right 
and Appendix A, while m^o starts 
yUc, then increases 



very slowly between and approximately = /i ~ Afj/v2, and increases very quickly 
afterwards. El 

More quantitatively, we find that there are two different values of Uc, that we dub 
/ic that are uniquely determined in terms of Vdi Vt and the A's. When one of these 
two critical values corresponds to m^o^ the other will correspond to m™J^, and vice versa, 
depending on the sign of the following quantity: 



fJ'C 



4(4Ai ~ A2 - A3 + 4X4)vf)/V'2{vj + 16vf). In fact ^ is the common 
value of /i at which the slopes of the variations of rn^o and rn^o as functions of /i experience a sudden 



^ The precise value is p, — vt{Xv'^ 
e c 
change 
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Va = (-A + Ai + A4) vl + A (A2 + A3) v''^ (7.1) 

Moreover, it turns out that at these extrema one of the two hP or states will correspond 
to a purely SM-like Higgs state, and this too is controlled by the sign of V\. One can 
summarize the behavior analytically as follows: 

i) Va>0 : 

In this case m/jo reaches a maximum given by 

2 max 2 _ tiy 

^ho (1) = "^(1) = ^ (7-2) 

when II takes the value 

// = /i«^(Ai + A4)^ (7.3) 
and m^o reaches a minimum given by 

= ^ + 0{v^) (7.4) 
when 11 takes the value 

^^^c'^^ /^,2'!i« ((2A-Ai-A4)^^ + 8(2Ai + 2A4-A2-A3)^,^). 
V2 (^2 + 16^2) 



(7.5) 



Expanding the Higgs masses squared around /j,^^ one finds 



mlo = m^o-- - + 0(5^,) (7.6) 



2 2 



Al,o = (-A + Ai + A4)^ + 2 (A2 + A3) + -j^5,, + 0{5l,) (7.7) 

A^o = (-A + Ai + A4)^ + 2 (Ai + A4) vi + y:^L_^s^^ + 0(52^) (7.8) 

= (-A + A, + ^)| + (2A, + A4) I + ^^^^^ V + 0{5l^ (7.9) 

Ah±± = (-A + Ai)^ - A3.,^ + -^5,, + 0{5%) (7.10) 



where = m\ ~ '^'hP^^ denotes the various squared mass sphttings from m^o™^ 
and 5fj_i = n — ^liP . 
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ii) Va<0 : 



In this case the reversed configuration occurs, m/jo reaches a maximum, given by 



2 max 



"^/lO (2) 



m 



(2) 



(7.11) 



(2) 

at = /ic , while tjiho reaches a minimum given by 



2 mill _ 2 
"^HO (2) ~ "^(1) 



(7.12) 



at /i = /ii^'', where m^^j^ym^^^y are as defined in Eqs.(|7.2 

Again, expanding around jjL'p we find 



7.5) 



2 max , ^"^ff [-2 , ,'/-\/ r-i \ 



(7.13) 



and the squared mass sphttings 



A^o = (A - Ai - A4) ^ - 2(A2 + A3) vl + -^5,2 + 0(5; 

2 V2wt 



(7.14) 



AO 



( (A - Ai - A4) ^ + 2 (2 (A - A2 - A3) - Ai - A4) vl ) 
+ ^ V + C'( V' /^rf) 



^"^ + '"'U,2 + 0(5j2,-f/-J) 



A 2 ^ 



((A - Ai - 2 A4) ^ - (4 A2 + 5 A3 + 8 A4) vl 



2 

-'1,2 



(7.15) 



11 



((A - Ai - - A4) ^ + (2 A - Ai - 4A2 - 4A3 - — A4) vf) 



(7.16) 



(7.17) 



where 5^2 = A* — 



(2) 



Noting that fi^c^ — \ m^2) ~ ^"^^^ ^\ have the same sign, and defining 



/zr^^max{/i«,/x(2)} 



(7.18) 
(7.19) 



one can recast the results of Eqs.(7.2, 7.4, 7.11, 7.12) in a more compact form as 
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mlo""^ = mlo {fx = /ir") = min{mji) , mf^) } (7.20) 

m^o™" = mjjoifJ, = /i™") = max{m^i),m(2)} (7-21) 

with an implicit reference to the two regimes i) and ii) if one keeps in mind that m^^^ is 
reached for /x = Hc^. 

The mixing pattern: for fi = fii\ hP and become pure doublet or triplet states, since 



in this case i? = as can be seen from Eq.(2.22). However, a close inspection of Eq.(2.40) 



shows that in the regime i) (resp. ii)) one has Sa = (resp. Sa = 1) for this value of /x. 
Thus, at /i = /ic^\ becomes a pure SM-like Higgs in regime i), but it is that becomes 
a pure SM-like Higgs in regime ii). The fact that the SM-like state is not always associated 
with the lightest CVeven state is important when discussing the Higgs phenomenology and 
interpretation of the experimental limits and is consistent with the fact that m^j^-j is indeed 



the SM Higgs squared mass, Eq.(7.2). In fact, due to the smallness of Vt/vd the behavior 
of the mixing angle a over the full range of the /i parameter follows closely the generic 
pattern discussed above: in both regimes i) and ii) one has essentially Sa — ±1 or Sq, ~ 
over most of the fx range, except for a very narrow region in the vicinity of /i defined in 
footnote [9] and satisfying 

n = \{^^^^'^ + ^^?^) (7.22) 

where |sq,| changes quickly from ~ to ~ 1. The generic dominance of no-mixing 
regimes can be understood from the asymptotic behavior at small and large /i values, 
i.e. I sinai^^ol = 1 - 2^-^^^^{vl/vl) + 0{v^t) and sina,^^^^ = 2{vt/vd) + 0{vl), together 
with the fact that dsa/d/x = 0{vf). We illustrated this behavior in Fig. [5] adopting the 
sign convention ea = +1. As seen in Fig. [5jb, remains positive in all the /i range since 



B < (cf. Eq.(2.22, 2.40)). And in accordance with the asymptotic behavior, Sa tends to 
(9(10^^) at large /i(> Jl) where is nearly SM-like, and to C(l) at small /i(< fi) where 

is nearly SM-like. [Note that in this numerical example /ii^"* becomes negative and is 
never reached]. In contrast, for the regime illustrated in Fig. |5}a, Sa remains negative for 
fi < fi^c \ crosses zero at fx^c^ and again tends to a positive value C(10~^) for ix ix^\ 

The exact magnitude of |sq,| at the three critical values of /i can be summarized as 
follows: 
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Va>0 



Va<0 



,(/x = 4^))|= , 1 

\ 1/2 / \ 1/2 

2 



(7.23) 



d 



(2)^1 ^ ^'d 4t't 

Large mixing scenarios have been discussed previously in [HI |38] while here we quantify 
more precisely the regions where such a large mixing takes place. 

For later analyses it is useful to characterize the n range in the large \sa\ regime. One sees 
from the above equations that the size of this range is 0{vt). As a first approximation one 



can characterize it by the interval < yU < /i™", with yU™™ given by Eq. (7.18). However, 
depending on the values of A and Ai + A4, \sa\ can still be very close to 1 in the range 
/i™" < fi < fl, especially that is not positive definite [it becomes negative when 
Ai + A4 < or 2A — (Ai + A4) < 0]. It is more sensible to base this characterization on the 

amount of deviation from the value \So,\ = 1. Defining fi in the vicinity of /2 in the form 

^2 ^3 

fi = p, — 6xVt, with 6 strictly > 0, one finds |sQ,(/i)| = 1 — k(6)-^ + O(^). For each given 

^d 

positive value of k there corresponds a value of fi given by 



/A A — Ai — A4 , s /_ .X 

The two-fold ambiguity in this expression is resolved as follows: requiring consistently 
^min < < ^ to hold, one should take for V\ > 0, fi = fi^-) with > 8, and for Va < 0, 



= fl[+) for any /c > In particular /t reproduces respectively /li^'* and for the 



special values k = and A; = 8 as expected, while fl cannot be reached for any finite value 
of k [consistently with the fact that |sQ,(/i)| ^ 1 and |sa(/i)| ~ I/V2 are not perturbatively 
close to each other in terms of powers of Vt/vd\- 

With the above prescription one can characterize the range in the large \sa\ regime 
by < /i < fi{k), where k can now be interpreted as triggering the experimental sensitivity 



to the deviation of \sa\ from its maximal value \sa\ = 1. Equation (7.24) shows that the 
lower the sensitivity to large \sa\ (i.e. the larger k), the lower the sensitivity of the size of 
the II domain to Ai + A4. We will come back to the above issues in the phenomenological 



""^^ Strictly speaking, in the case Va > one can still choose k in the interval 2 < /c < 8 if Va is sufficiently 
close to zero so that to ensure that < /t. In practice these details will not be important, since one 
does not expect an experimental sensitivity to the deviation from |sq| = 1 to be better than a few percent. 
A deviation of 1%, with vt = IGeV, puts the value of k already around 600 ! 
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Figure 3: Illustration of the regime Va < with A 
-i,A4 = l,vt = l GeV,v = 246 GeV,^;^ = ^yv^^-2vlK 
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IGtt \ 



2' 



10-\A3 = 2 X 10-\Ai = 
8, leading to fi'c^ ~ 23 GeV. 



see Eq. (7.5). 



discussion of section |8l 

Unitarity bounds: relying on the above properties we can now easily derive the theoretical 



upper bounds on the various Higgs masses. From Eq. (7.2), and using the maximal value 
allowed by the tree-level unitarity constraint for A, Eq. (6.14), and Vd — 246GeV 
determine an upper bound on m/jo. 



we 



niho < 712 GeV (for k = 8) (7.25) 
< 1 TeV (for k = 16) (7.26) 



If Eq. (7.11) is used instead, then the saturation of unitarity and BFB bounds on Ai + A4 
should be also considered. However, due to the smallness of Vt/vd, this would lead to 
only a few GeV change in the above upper bounds. As far as m^^o is concerned, the 



above bounds are essentially the minimally allowed values, as obvious from Eq. (7.20 



7.21), in the unitarity saturation limit. To obtain its theoretical upper bound as well as 
those of the other Higgs masses, one should rather take n at its maximally allowed value, 
— /i+ — -^^i;^} since all these masses increase monotonically with /x. For instance, 
with the set of parameters chosen in Figs. |3l |4] and Vt = I GeV one finds the upper bounds 



mH±± — rnH± ~ — rriHO ~ 88 TeV (for k = 8) (7.27) 

~ 124 TeV (for k = 16) (7.28) 

which are not phenomenologically compelling. Actually, somewhat lower bounds are ob- 
tained when taking into account the experimental exclusion limits on a light Higgs m/jO but 
remain too high to be useful. In contrast, phenomenologically interesting scenarios with 
light charged, doubly charged, CVodA and CVeven Higgses are possible for small values of /i. 
For instance as illustrated in Figs. |3| |4| such a light spectrum occurs when n <^ /xi^^ ~ 23 
GeV. More generally, the analytical expressions given above for the mass splittings show 
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that in the vicinity of /^c and in particular for [i < [i^ the neutral CP even is not necessarily 
the hghtest Higgsp] The detailed patterns will depend on the actual values of the A's and 
will be studied more thoroughly in the next section, but one can already see some generic 
features in regimes i) and ii) at /x ~ /ic. In regime i) where —A + Ai + A4 > one expects 
if^^ to become the lightest Higgs if — A + Ai < 0, that is when Ai < A < A1 + A4. Similarly, 
in regime ii) where typically A — Ai — A4 > one again expects H"^"^ to be the lightest 



Higgs when Ai + A4 < A < Ai + 2A4. More generally, a close inspection of Eqs.(2.11, 2.27^ 
shows that mH±± < m/jO when /i < /i* = (A + X4)vt/\/2 + 0{v^/v'^), and only if A4 > 



Furthermore, it immediately follows from Eqs.(2.13, 2.31) that mH±± < rn^i < rnA when 
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A4 > so that the necessary and sufficient condition for H^^ to be the lightest Higgs is 

fi< fi* with A4 > 0. (7.29) 
Phenomenological bounds: in order to prepare for a phenomenological study, we discussed 



in section 3.2 the modification on the tachyonic bounds of fi when experimental exclusion 
limits are available for m^o,m/^± and mH±±, cf. Eqs.(3.65, 3.66). Here we address the 



same question concerning m/jo and 111^0 ■ Given some experimental exclusion lower bounds 
(m/io)exp, (resp. (m/^o)exp) there corresponds two values (resp. /i±°), namely 



8V2vt 
± 2 [(mji) 



"^ho)exp)("^?2) - 



cxp 



hO Jcxp, 



)]H1 + 16 



t \i 



(7.30) 



(and similarly for fi^° with (m/jo)cxp replaced by (m//o)exp ), for which m/jo reaches (m/jo)exp 
(resp. rriHO reaches (wi^'Oexp)- Note that in the limit of no experimental bounds, i.e. 



^^We have kept in these expressions subleading terms of 0{v'^) in order to handle as weU the small parts 
of the Ai's parameter space where the leading 0{v'^) are suppressed. 



^'^ Although this expression of /i* is well-defined for A4 < 0, one finds that the splitting m 
is negative only in the domain defined by /i < (A + 2X4^)vt/\/2 + 0{vf /v"^) and — A4U^/(4A/2wt 
/I < (A + X4)vt/V2 + 0{v'l/v^), which is clearly non-empty only for A4 > 0. 



2 

_f/±± 



Oivt) < 
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— 0, Eq.(7.30) gives back Eq.(A.l). Furthermore, relying on the fact that 



m^o has a maximum and m^^o has a minimum as functions of /x, cf. Eq. (7.20, 7.21), the 
phenomenological bounds read 



IJ''^' < < /^+ assuming (m/io)exp < m 
and 



max 
hO 



(7.31) 



/i < yu;'^'' or < yU assuming {rnHo)exp > w^^o" 



Obviously (m/io)exp > m^^^ would be an inconsistent assumption, while (mi^o)exp < rn 



''HO 



would be an empty assumption not leading to any constraint as far as /i is concerned. 
In summary, the experimental lower bounds on the various Higgs masses will typically 



constrain the /i parameter to lie in a finite domain defined by the combination of Eqs.(3.66 
7!3T|. 



8 Higgs phenomenology 

Although previous studies in the literature assumed typically the triplet mass Ma and the 
mass parameter fi to be much larger than the electroweak scale, ^ Vd, attention has 
been paid more recently to the possibility of having M/^, fj, < 1 TeV where the Higgses of 
the DTHM might be accessible at the Tevatron and the LHC [HI EH SHI SD 1121 SSI IH 
H5l HE]. In this spirit, the results obtained in the previous sections help defining educated 
strategies to extract constraints on the physical Higgs masses and model parameters from 
experimental data, rather than performing merely blind (and CPU time consuming) scans 
on these parameters. The existing experimental exclusion limits on the SM Higgs particle 
are readily translated into constraints on the DTHM in the parameter space region where 
becomes SM-like, i.e. when the mixing between the doublet and the triplet is very small. 
However, even when far from this region, existing exclusion limits for an extended Higgs 
sector (such as in two Higgs doublets models or in the minimal supersymmetric extension 
of the SM) can also be partially adapted to h^, , and H^, while of course H^^ has 
a distinctive experimental search. 

In this section we give a quick overview of the Higgs sector phenomenology and exper- 
imental searches (for an extended overview on the phenomenology of triplet models see 
Ref. jlZ]), followed by a preliminary analysis using our results. A detailed study taking 
into account all present-day experimental limits lies out of the scope of this paper and will 
be presented elsewhere. 

Doubly charged Higgs 

Observation of the doubly charged Higgs H^^ would signal unambiguously physics beyond 
Higgs doublets, let alone physics beyond the SM Higgs sector. Owing to charge conserva- 
tion, it is obvious that H^^ cannot couple to a pair of quarks, therefore, its possible decay 
modes are: 

i) same sign charged lepton pair H^^ — )■ l^l^ that proceed via lepton number violating 
coupling, 

ii) a pair of gauge bosons H^^ — )■ 

iii) H^^ W^H^, 
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iv) a pair of charged Higgs bosons H^^ — )■ H^H^ . 



We emphasize also that the doubly charged Higgs couples to the photon and to the Z 
boson through gauge couplings Eqs. (C.17, C.18), while its couplings to a pair of is 
proportional to the triplet vev fj, see Eqs. (C.15). Therefore the decay channel ii) will be 
suppressed for vt ^ Vd- The decay channel iv) will also be suppressed for small Vt as can 
be seen from the form of the coupling of H^^ to a p air ch arged Higgses H^, Eq. (C.16). 
Indeed, one has cos /3' ~ 1 and sin /3' ~ Vt/vd from Eq. (2.19), and furthermore, the /xsin 
is also of order vt due to the /i upper bound /i+ ~ Vd/'^'t, viz Eq. (A.l). In contrast, the 
coupling H'^'^W'^H'^ is proportional to the gauge coupling and has no suppression factors. 
The decay channel ii) will thus contribute substantially if kinematically open. Depending 
on the size of the Yukawa couplings of the leptons, the doubly charged Higgs can decay 
dominantly either to a pair of leptons or to and H^, and subdominantly to a pair of 
and/or a pair of if kinematically allowed. 
In e~^e~ collisions, the doubly charged Higgs can be pair produced through 7 and Z 
s-channe|^ e+e" 7*, Z* H^^H"^^ [HI HHl |50l EI] . One can have also access to the 
associate production of with through s-channel Z exchange [52t l53t [5l] . If the 
e~e~ option is available at ILC, then the doubly charged Higgs can be produced in W^W^ 
fusion through e~e~ — )■ W~*W~* — )■ e~e~H^^ . Even if H^^W^W^ has a Vt suppression, 
the rate for W^W^ fusion could be substantial especially at higher energies options for 
e-e" [19]. 



At the Tevatron or the LHC, the two production mechanisms with potentially large 
cross sections are pp/pp — )■ 7*,Z* — >■ H^^H^'^X or a single production through WW 
fusion pp/pp — )■ — )■ H^^X [551 ES]. The latter process as well as the s-channel 
pp/pp — )■ W^* — )■ W^H^^ depend on the coupling H^^W^W^ which is proportional 
to the triplet vev. However, the suppression due to the small value of Vt is somewhat 
compensated by the fact that IV^IV^ fusion could be substantial at high energy. Those 
processes have to be supplemented by the associate production of singly and doubly charged 
Higgs bosons pp/pp — )■ H^^H^X which could have a comparable cross section to pp/pp — )• 
H^^H^^X [HZ], [58]. 

Such doubly charged Higgs have been subject to many experimental searches. At 
LEP-II, the experiments L3, OPAL and Delphi [59], [60], [61] performed a search for 
doubly charged Higgs boson assuming that H^^ decay dominantly to a pair of leptons 
H^^ — )• l^l^. Four leptons final states have been analyzed at L3, OPAL and Delphi. L3 
performed a search for the six possibilities: ee, /i/x, fie, fir, er and rr. No excess has been 
found and lower limits in the range 95-100 GeV at 95% confidence level on the doubly 
charged Higgs boson mass are derived. Those lower limits depend on the doubly charged 
Higgs decay modes. For example, if H^^ — )■ e^e^ is the dominant decay, then the lower 
limit is 100 GeV while if H^^ — )■ fi^r^ is the dominant decay then the lower limit is about 
95 GeV. 

At the Tevatron, D0[62],[63] and CDF [M], [65] have searched for ^ 7*, Z* ^ H^^H^^X 
with H^^ — )■ l^l^. D0 measurement [62] represents the first doubly charged Higgs search 
with the decay H^^ — )■ /z^/i^. Note that D0 search was limited to H^^ — )• /x^/z^ that is 
an almost background free signal, while CDF explored the three final states e^e^,/i^/i^ 
and e^/i^. Both D0 and CDF excluded a doubly charged Higgs with a mass in the range 
100 — 7- 150 GeV. We stress that all those bounds assume a 100% branching ratio for 
H^^ — )■ l^l^ decay, while in realistic cases one can easily find scenarios where H^^ — )■ l^l^ 

^^the t-channel mediated by a lepton is in general suppressed by the small Yukawa coupling 
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is suppressed while H^^ — )■ is substantial [171 |39l |66l [67] which could invalidate 

partially the CDF and D0 limits. However, the LHC has the capability to extend the above 
limits up to a mass about 1 Te V for high luminosity option [TTl HOj |56l [68] . Observation of 
doubly charged Higgs bosons at the LHC and measurement of its leptonic branching ratios 
will shed also some light on the neutrino mass pattern [T71 [3S1 EH ISSl [SSI [SZl [SHI [ZOl [B] • 
Finally, indirect limits on the mass and the bileptonic couplings of the doubly charged 
Higgs boson can be extracted from low energy lepton flavor violating processes, such as 
/i — 7- 67, yU — 7- 3e, r — 7- 3Z,... (see for instance jH] 1^). 



Singly charged Higgs 

Let us now discuss briefly the couplings of the singly charged Higgs and its decay modes. 
The charged Higgs coupling to lepton and neutrino is proportional to ~ m^/vt ~ [T7] 
which could be of the order 0{1) if Vt is very small. Similarly, the charged Higgs coupling 
to a pair of quark u and d is proportional to tan/3' which is suppressed by Vt/vd [IZ]. In 
the case of H~ib, this coupling could enjoy some enhancement from Yukawa coupling of 
the top quark. The suppression of the coupling H~ib has three consequences: 

• Given the suppression factor of the order Vt/vd for H~ib, the charged Higgs mass can 
not be subject to 6 — )■ S7 constraint, similarly to the two Higgs doublet model type 
I where the coupling is suppressed by 1/ tan/3. 

• Some of the conventional mechanisms for charged Higgs production at Hadron col- 
liders such as bg — tH~^ and gg — > tbH^ will be suppressed. 

• Since the charged Higgs search at the Tevatron is based on the top decay t H^b, 
given the suppression of H~^tb coupling the branching ratio of t — )■ bH~^ would also 
be suppressed. One concludes then that the CDF limit does not apply in this case. 



Besides those processes which are suppressed, one can still produce charged Higgses 
through the associate production of singly and doubly charged Higgs pp/pp — )■ W* — )• 
H^^H"^ [571 EH] with a spectacular signature from H^^ — )■ l^l^. Other mechanisms are: 
the Drell-Yan process pp / pp ^ '-)* , Z* — )■ H^H^, the associate production of charged Higgs 
and neutral Higgs pp/pp — )■ W* — )■ H^h^, pp/pp — )■ W* — )• H^H^, pp/pp ^ W* — )■ H^A^ 
and the associate production of charged Higgs with W gauge boson pp/pp — )• Z* — > W^H^. 
Note that among the latter processes, the ones with W^H"^ or H^h^ final states are 
suppressed by a Vt/vd factor as compared to the Drell-Yan and the two other associate 



production processes that are controlled by gauge couplings, cf. Eqs. (C.IO, C.13) and 
Eqs. ( [aITj[Cl2| . 

If the charged Higgs decays dominantly to leptons (for small Vt) we can apply the LEP 
mass lower bounds that are of the order of 80 GeV [7T], [72]. For large Vt, i.e. much larger 
than the neutrino masses but still well below the electroweak scale, the dominant decay is 
either if+ tb or one of the bosonic decays if+ W~^Z, W^h^ /W^A^ . For the 

first two decay modes there has been no explicit search neither at LEP nor at the Tevatron, 
while for the — )■ decay (and possibly for — )■ W^h^ / if decays similarly 

to one can use the LEPII search performed in the framework of two Higgs doublet 
models. In this case the charged Higgs mass limit is again of the order of 80 GeV [72] . 
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Neutral Higgses 



The lighter CVeven Higgs boson is fully dominated by the doublet component (i.e the 
mixing \sa\ <^ 1) when fx > fx, as discussed in section [T] and illustrated on Fig. |5} In this 
case the coupling of to a pair of neutrinos is suppressed being proportional to Sa- Such 
Higgs will completely mimic the SM Higgs boson and then the LEP and the recent Tevatron 
limits would apply. In this scenario of very small mixing, the other neutral Higgses and 
would be fermiophobic to all charged leptons and quarks but their coupling to a pair 
of neutrinos that is proportional to cos aY^, k-Yj^ = niy/vt could be enhanced for small Vt- 
Then the dominant decay mode for and ^4°, for small Vt-, would be a pair of neutrinos 

Note that A^, being CVodd, does not couple to a pair of gauge bosons while the couplings 
H^ZZ and H^WW in the small mixing case are suppressed by Vt/vd-, Eqs.(C.4, C.6). 
Thus the W and Z Higgsstrahlung productions of and A^ are expected to be small. 
Furthermore, while the H^A^Z vertex is controlled by the gauge coupling, Eq. (C.8), h^A^Z 
has an extra Vt/vd suppression, Eq. (C.7). This implies that in the small mixing case, one 
can still produce A^ and through the Drell-Yan process e^e~ /pp/pp — t- Z* — t- H^A^. 
For very small Vt-, A and would decay essentially into a pair of neutrinos. At LEP, 
the signal would then be a photon (from initial state radiation) and missing energy in the 
final state. A lower bound on uih and of the order of 55 GeV can be extracted in this 
case from LEPII data, assuming mass degeneracy between A° and [73J. [In the non 
degenerate case the lower bound translates into mu + > 110 GeV.] Increasing vt well 
above the neutrino masses decreases significantly /A^ — )• uv, and the decay channels 
hi, A^ hh, as well as ZZ, A^ Zh^ if open, see Eqs. ( |C.4[|C.7[ ), can become 

dominant. If — bb, A^ — )■ bb dominate, the LEPII Higgs search through e^e~ — ?• H^A^ 
in the two Higgs doublet Model can apply to the DTHM in this case, and the limit is 
roughly tuh + TnA > 185 GeV [71]. There is however a distinctive feature in the DTHM 
related to the H^h^h^ coupling, Eq. (C.9), the latter becoming substantial for increasing 
/i and thus for heavier H^. The — )■ h^h^ decay mode would then be important in both 
the large and small (with respect to the neutrino masses) Vt regimes. 

In the case of maximal mixing \sa\ ~ 1 which occurs for jj, < Ji, see Fig. [sj the roles of 
and if ° are interchanged. is fully doublet and is fully triplet. Taking into account 
this interchange, the previous discussion applies here to . However, since remains 
the lighter Higgs which can be now far from SM-like, one expects weaker experimental 
constraints on its mass than the ones quoted above. 



Top decay into charged Higgs 

A light charged Higgs of the order 100 — 200 GeV is still allowed by theoretical constraints 
as well as by experimental search. If the charged Higgs satisfies mH± < mt — rrib, one could 
ask whether the decay t — )■ bH^ can have a significant branching ratio to be observed at the 
LHC. As mentioned before, the coupling H^tb has a Vt/vd suppression and the branching 
ratio for t — ?■ bH^ is expected to be small. We perform a systematic scan over the DTHM 
parameters space looking for charged Higgs masses that allow the t — )• bH^ decay to be 
open. In Fig. |6] we show the branching ratio for t — )• bH^ where we included t — )• bW^ 
and t — )■ sW^ decay channels and the QCD corrections. It is obvious that a large effect 
on t — )• bH^ would appear for the largest possible values of Vt that are allowed by EW 
precision constraints and the theoretical constraints. Indeed for a triplet vev vt in the range 
0.1 — )■ 3.5 GeV and a charged Higgs mass less than 165 GeV, one finds Br(t — )■ bH^) in 
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II (GeV) II (GeV) 

Figure 5: The mixing angle as a function of fi, in the regimes Va > (a) and Va < (b); 
the other parameters are given by t>t = 1 GeV, Xmax = 167r/3, A2 = A3 = 0.1, Ai = 0.5 and 
ea = +. The log scale in (b) shows the asymptotic values at large fi. The same asymptotic 
values apply in (a); see text for further discussion. 

the range 10"^ 10~^. 

However, it is well known that the LHC will act as a top factory. With low luminosity 
10fb"\ 8 million ti pairs per experiment per year will be produced. This number will 
increase by one order of magnitude with the high luminosity option. Therefore, the prop- 
erties of top quarks can be examined with significant precision at LHC. For instance, it has 
been shown that for top decays through flavor changing neutral processes, it is possible 
to reach Br(t — )• cif°) < 4.5 x 10~^ at the LHC [75]. For t — )■ bH~^, no such studies are 
available. But it is clear that if we let one top decay to bW and the other one decay to 
bH with a branching ratio in the range 10^^ — )■ 10"'^, this would lead to 800 — )■ 8000 raw 
bW~^bH~ (or bW^bH^ events in the case of high luminosity option which may be enough 
to extract charged Higgs and measure its coupling to the top. Note also that high sensitiv- 
ity to the charged or neutral Higgses of top decays through loop induced flavor changing 
neutral currents, can also be attained at the ILC [761 EH EH] . 



DTHM spectrum and theoretical constraints 

We illustrate, in Figs. [7}a and|8}a, the correlations among /i, sin a and Vt for fixed values of 
the Aj's and A, and in Figs. [7}b, [8jb and[8jc, the correlations among /i, sin a and the CVeven 
Higgs masses (or equivalently A), for fixed values of the Aj's and Vt, where we take into 
account the boundedness from below and unitarity constraints discussed in the previous 
sections. Note that the chosen numbers in the figures are such that Va < in Figs. [7| 
and Va > in FigMa, while FigsjSjb and [8jc interpolate between these two regimes, see 



Eq. (7.1). For fixed /x, increasing the magnitude of Vt decreases ruho and increases the 



mixing parameter |sq,| as can be seen from Figs. [7|a andjsja. The upper-left white areas 
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Figure 6: Branching ratio for t — t- bH^ 
values of the triplet vev Vt- 



as a function of charged Higgs mass for three 



in these plots correspond to m/jO < llSGeV where we took the latter value as a fiducial 
lower bound for a standard-model like Higgs. Such a bound corresponds to A^'^ ~ 0.44, 
for Vt < IGeV, while the upper bound for ruho is around 120GeV, corresponding to the 
value A = 0.48 chosen in the figures, cf. Eqs. (7.2, 7.4). It thus follows that the colored 



areas in the plots, indicating mainly very small values, i.e. behaving like a SM Higgs, 
correspond to the small Higgs mass range llSGeV < m/jO < 120GeV. Increasing the value 
of A, keeping A^^ fixed, would result in an increase of the Higgs mass range as well as 
of the regions with larger \sa\ (the red areas on the plots). In fact, there are two regions 
corresponding to rriho < llSGeV, the white area in the upper-left corner corresponding to 
small values of fi delimited by the red thin area, and another region at very large values of 
fi (> (9(1) — (9(10'^)TeV), which are out of the scope of the fi range shown on the plots, that 
are delimited by green-blue areas. One should note that, in the former region, \sa\ reaches 
quickly 1, so that carries essentially a triplet component and is thus not excluded by a 
fiducial SM-like Higgs mass lower bound, even if it is lighter than this bound. In contrast, 
in the latter region where \sa \ remains very small, a SM Higgs mass lower bound applies to 

It follows that such a bound does not put lower bounds on /x, while it leads typically to 
very large upper bounds on as a function of Vt. In the small /i region, carries mainly 
the SM-like component and should respect a SM Higgs mass lower bound. However, due to 
the very low sensitivity to fi in this regime (see Fig. |3]), such a bound will translate merely 
into a lower bound on A. Therefore, exclusion of very small values of fi can only originate 
from exclusion limits on the lightest non-SM-like CVeven or CVodd Higgses, which could be 
extracted for instance from existing limits for the the minimal supersymmetric extension 
of SM in the non-decoupling regime [79] . 

Complementary features, now with a fixed vt and varying A, are illustrated in figures 
[7|b, [8jb and|8|c. The gross features of Figs. [7|b and|8}b are in agreement with the previous 
discussion on the phenomenological bounds, related to Eqs. (7.30, 7.31). They illustrate 



how an information on m/jO constrains the allowed range for /i without any prior knowledge 
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on A. For a given m/^o, the allowed range of A is theoretically bounded from below by some 
Amin, in order to satisfy m^o < m™o^^, see Eq. (7.31). Then for each value of A in the 



domain Amin < A < Amax = there corresponds two values of /i consistent with a given 



rriho, according to Eq.(7.30). Then it is easy to see, from the shape of the m/io(/i) plots 
shown in Fig. pi that the largest spread between and fi'^' is reached for A = Amax, since 
increasing A results in shifting upwards these plots . The two branches of the envelop of 
the domains in Figs. [Tjb, |8|b correspond to fi'f (Xmax) ■ Furthermore, increasing m^jo with 
fixed A = Amax results in an increase of /i_ and decrease of as can be again seen from 
the shape of mho{n) plots shown in F ig. [3 till the two branches join and terminate when 
rrifiO reaches its unitarity bound Eq. (7.26). With the numbers chosen on the plots, /i is 
bounded to lie between /i_ 0.3 GeV and /i+ ~ 10^ GeV. One can see that for small /i < 1 
GeV, rriho must be less than about 200 GeV. The latter bound on rriho increases quickly 



to reach the unitarity bound Eq. (7.26) when /i increases from 1 GeV to 10 GeV. Above, 
/i = 10 GeV, rriho can be any number between the LEP limit (114 GeV) and this unitarity 
bound. As noted previously, one should take into account the actual doublet content of 
when reading out exclusion domains from these plots. In the plot, we have illustrated the 
size of |sq,|. In most of the cases the mixing angle is very small (blue to green areas), which 
means that is dominated by doublet component. In these regions where a SM Higgs 
exclusion limit can be readily applied, one might still need to combine this information 
with the search limits for the other charged, doubly-charged and CVodd Higgs states, in 



order to reduce further the otherwise large allowed domain for /x, see Eq. (3.65). However 



due to the Vt suppression in Eq. (3.65) of the lower bound /imin, such a reduction is not 
expected to be significant unless the experimental lower bounds, (m//±±)cxp or (m//±)exp 
or (m^Ao)exp, become sufficiently higher than the electroweak scale. In contrast, bounds 
on m/jO alone would shrink significantly the spread of the /z range whenever \sa\ > 10~^ 
(the green/red areas), reducing as well the order of magnitude of the size of fi. In such a 
regime of small n one starts being sensitive to values of the Aj's, as can be seen through 
the slight difference, in the green area, between Fig. [7|b and Fig. [Sjb. This effect will of 
course increase for higher values of the Aj's consistent with unitarity and BFB constraints. 

Figure [8jc illustrates the behavior of rriHO as a function of fi and A which, as compared 
to Fig. [8jb, shows a striking difference from the behavior of m/jO. According to the pre- 
vious discussion on neutral Higgses (see also Fig. [s]), \sa\ is essentially either very small 
or very close to 1. Thus, the red area corresponds to an behaving essentially like the 
SM Higgs. The dual sizes of the red areas in both plots can be understood again from the 
mass shapes of Fig. [3] for small /i(< fx), m/jO changes very quickly with /i while 711^0 is 
almost insensitive to /z. It follows that a variation of A, that amounts to shifting upwards 
or downwards these mass shapes in Fig. [3} results in a small change in /i for a fixed m/jO 
and a big change in fi for a fixed itlho, whence the narrow red strip in Fig. [Sjb and the 
large red area in Fig. |8}c. (One can understand similarly the dual sizes of the blue and 
green areas for large /i(> p,).) These features suggest a useful complementary strategy 
when using present or future exclusion limits, depending on whether one interprets these 
limits in the small or large \sa\ regimes. We discuss this strategy only qualitatively here, 
summarizing its main points as follows: 

-I- in the small \sa\ regime, akin to moderate to large /i values, the typical Higgs spectrum 
features a CVeven lightest state /i° behaving like a SM-Higgs, the remaining Higgs states 
being much heavier as illustrated in Fig. [31 and Fig. |9la. Interpreting the exclusion limits 

— — (SM) 

within this regime amounts to applying a SM Higgs mass lower bound m\ to m/^o that 



leads to a lower bound on A, see Eq. (7.31). To any A above this bound will correspond a 
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lower, , and an upper, /x^", bound on /i. The lower bound is, however, typically too 
small to be consistent with the small \sa\ regime and should be superseded by a larger value 
(!?(max{yUc^'*, /ic^ '')}). Furthermore, one should keep in mind that is extremely sensitive 
to rriifi and decreases quickly with increasing m/jo. This implies the important feature 
that a slight improvement of the exclusion limit mf^'' results in a substantial decrease of 
the upper bound on fi. The heavier CVeven state is not expected to bring significant 
constraints. Indeed, in the considered regime, this state carries essentially the triplet com- 
ponent with suppressed couplings to the SM sector. Its mass can thus be bounded only by 
^(non SM)^ exclusiou mass limit on non-SM-like Higgs particles. Since such an exclusion 
mass limit is expected to be weaker than the SM-like limit due to lower statistics, that is 
m|j^°"~^^^ < m!^^\ then taking into account that one has theoretically mfjo > rriyfi^ one 
is trivially lead to m^^^"^ ^^'^ < m^™ which implies no new constraints (cf. the discussion 



following Eq. (7.31)). As stated previously, exclusion limits on the remaining Higgs states 



can also be used independently to improve the lower bound on based on Eq. (3.65). One 
can, however, get further information within the present regime depending on whether 
these exclusion limits are higher or lower than m^^^\ In particular if {mH±±)exp ^ "^If^^ 
which excludes an H^^ lighter than h^, then one excludes all the A4 > region, or else. 



puts a stronger lower bound on /x, namely /i > /i*. (see Eq. (7.29) and discussion thereof.) 
In the case where (m//±±)cxp ^ ^^h^\ which is the present experimental situation, there 
is a small window firnin < fJ^ < fJ^* with A4 > 0, otherwise ^* < fi < fi'^ irrespective of 
the sign of A4, and for all the allowed values of A discussed above. We have illustrated in 
Fig. [ojb a case where H^^ can be the lightest Higgs state. 

-II- in the large \sa\ regime, akin to small fi values, is the heaviest among all the Higgs 
states of the model and behaves like a SM-Higgs. This is a rather unusual configuration 
that should help constrain more efficiently, or perhaps exclude, this regime. Also in this 
small /i regime, and in contrast with the previous regime where only A was playing a role, 
there can be now a somewhat increased sensitivity to the Aj's as well, in particular Ai + A4. 



The reason is that the size of the fi domain is of order ft, Eq. (7.24), where in the latter 
Ai,A4 do not suffer a Vt suppression as compared to A. However, as discussed in section 
7\ the parameter k will characterize the sensitivity to the deviation of from a pure 
SM-Higgs state, which can lead, for realistic experimental sensitivities, to a significant 
reduction of the sensitivity on Ai + A4. 
One then has to consider two cases: 



a) m^h^^ < m™o°: this case implies essentially a lower bound on A through Eq. (7.21 ), but 
no constraint on /z apart from the defining region in this regime, namely < fi < fi, 
whose size depends mainly on A and to a lesser extent on Ai -t- A4. The latter couplings 
are bounded by the combined unitarity and BFB constraints of section [6} so that 
there is an indirect sensitivity to A2 and A3 as well. The red area in Fig. [Sjc gives 
an illustration of this least constrained case. The fi domain extends over all the red 
area, while the vertical boundary of this area is determined by the maximal value of 
A = ^ given by unitarity. This boundary corresponds to the unitarity upper bound 



on the SM-Higgs mass as well as the one on m/jO, Eq. (7.26). Of course can escape 
this bound but at the expense of switching consistently to the small \sa\ regime as 
seen on Fig. [8jc. 

t*) "^/f^'' ~ m™^^: in this case not only do we have an upper bound on A through 



Eq. (7.21), but actually also a lower bound. Indeed, a too small A, leading to a 
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too low m^o" with respect to rn^h^\ will eventually put rn^^^^ J^^^ above all the 
values of rrifjo corresponding to the large \sa\ regime thus ruling out this regime 
altogether. Furthermore, this configuration will immediately rule out the small \s 



a 



regime as well, since m^^^^ is by definition applicable only to SM-like states and we 
have m^h^^ > m^™ > ruho for all m/jo. Consequently, a too small A would exclude 
the whole fi parameter space. One concludes that A should lie in a very narrow 
strip such that m^™ < m^h^^ < rriHoin = 0). This strip is essentially giving the 
lower bound on A of case a) and thus does not provide significantly new information. 
[Note, though, that for values of Ai + A4 close to its unitarity bound, and taking for 
instance m^h^^ — 114GeV and vt = IGeV, case b) can still reduce the lower bound 
on A, from A ~ 0.43 to A ~ 0.38. But the effect will be smaller for smaller values 
of Vt] One should, however, keep in mind that due to the high flatness of m^/o as a 
function of /i in the region < /i < /(, the slightest variation of A within the above 
noted strip would result in the exclusion of signiflcant parts of the < ^ < fi region. 
For instance, with Ai = 1.5, A2 = A3 = 0.1, A4 = —1 and taking the SM-Higgs lower 
bound mf^^ = 114GeV, if one reduces the lower bound of A (for which the whole 
range < < /i ~ 0.3 of the large \sa\ regime is allowed) by just l%o, then the 
SM-Higgs lower bound would imply /i < — 5GeV or > 0.3GeV thus ruling out the 
whole large \sa\ regime! 

For all practical purposes and barring the flne-tuned effects just mentioned, the above 
discussion of cases a) and b) shows that an experimental lower bound on the SM-like Higgs 
mass cannot by itself cut into the large |so| /small fi regime; it either excludes it or allows 
all of it, depending on whether A is respectively below or above the value A that satisfles 
m{f^^ = m^jf(A = A). Thus, the size of the n domain [0,fi] that is controlled mainly 



by A [but can also be sensitive to Ai -t- A4] for each given value of Vt, see Eq. (7.24), will 



not be reduced by the actual value of the experimental bound m^^^\ Moreover, an extra 
constraint from an experimental lower bound on the mass of a non-SM-like CVeven Higgs 
state would have a marginal effect since m/jo decreases very quickly in the region ^ ^ fi- 
An efficient reduction of the /i domain can come only from experimental lower bounds on 
the masses of the charged, doubly-charged and CVodd Higgs states. Indeed, these bounds 



translate into a lower bound on /i typically of the same size as fi , Eq. (3.65). As far as 
these experimental bounds are of the same order, (m^o)exp — (^j^±)exp — (?Ti^±±)oxp, the 
relevant bound will be given by (resp. H^^) when A4 < (resp. A4 > 0). 

Comparing /imin and /i, one determines easily the necessary and sufficient conditions 
for which the large \sa\ regime would be excluded. They read: 

(m^o)exp > (A;A - 2(Ai + A4) - V2k\\ - Ai - ^41)^^^^^^ + 0{vl), for A4 < 0, and 
(m^±±)exp > {k{X - A4) - 2Ai - y/2k\\ - Ai - ^41)^^^^^^ + 0{vl), for A4 > 0, where we 



have taken into account the two-fold structure as discussed after Eq. (7.24). 



9 Conclusion 

We have carried out a detailed study of the renormalizable Higgs potential relevant to the 
type II seesaw model, keeping the full set of the seven free parameters of the potential. 
We determined analytically the unitarity constraints on the various scalar couplings and 
fully solved the all directions conditions for boundedness from below. These combined 
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Figure 7: (a) correlation between fi and Vt with niho > m 



(SM) 



llSGeV and A = 0.48, 



(b) correlation between yU and m^jo, scanning over A in the range 0.44 < A < 167r/3, with 
Vt = IGeV; color code: lO^^ < s„ < 1 (red), lO^^ < < iQ-i (green), 10" 
(blue) and Sa < 10^'^ (magenta). The other parameters are A^ — — A4 — . 
and K = 8. Va < for both figures. 



A2 = 



< 10-2 

A3 = 



theoretical constraints delineate efficiently the physically allowed regions of the parameter 
space and should be taken into account in phenomenological studies. We also examined the 
vacuum structure of the potential and determined general consistency constraints on the /x 
parameter, as well as theoretical upper (resp. lower) bounds on the lighter (resp. heavier) 
CVeven Higgs particlc mass that can constrain further the phenomenological analyses. We 
also identified two distinct regimes respectively for large and small /i. In the first regime 
the lightest Higgs particle is the h^, behaving as a SM-like Higgs, the remaining Higgses 
being typically too heavy to be of any phenomenlogical relevance. In the second regime, 
it is the heaviest Higgs that behaves as a SM-like Higgs, the lighter charged, doubly 
charged and CVodd states become accessible to the colliders, with possibly the H^^ being 
the lightest state, while the lighter CVeven decouples quickly from the SM sector. We also 
initiated the study of possible consequences from existing experimental exclusion limits. 

Although we did not commit to any underlying GUT assumptions, thus allowing /i to 
vary between a few GeV and possibly the GUT scale, we do retrieve, as a consequence of 
the (model-independent) dynamical constraints on yU, a seesaw-like behavior that leads to 
tiny Vt if /i is taken very large. 

Finally, the results of this study having been obtained at the tree-level, one should 
keep in mind that modifications due to quantum corrections to the effective potential can 
possibly be substantial in some cases. The inclusion of such corrections is, however, beyond 
the scope of the present paper given the non-trivial form of the constraints already at the 
tree- level. 
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500 1000 1500 2000 2500 

iJiH" (GeV) 

Figure 8: (a) correlation between and Vt with rriho > mfj^^ = llSGeV and A = 0.48, 
(b) correlation between fi and the light CVeven Higgs mass, (c) between /i and the heavy 
CVeven Higgs mass, scanning over A in the range 0.44 < A < IGvr/S, with Vt = IGeV; 
color code: 10^^ < < 1 (red), lO^^ < j^^j < iq-i (green), lO^^ < \sa\ < 10~^ (blue) 
and \sa\ < 10~^ (white bottom area in (a)). The other parameters are given by Ai = 1.5, 
A2 = A3 = 0.1, A4 = —1 and K = 8. Va > in (a), while in (b) and (c) Va changes sign 
with increasing Higgs masses. 
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Figure 9: Higgs boson masses as a function of /i with Vt = I GeV, A = Svr/S, Ai = 0.5, 
A2 = A3 = 0.1, for Va < 0, A4 = -1 (left) and Va > 0, A4 = 10 (right), we note that in 
both panels mno = m^o. 



Appendix A. 

As stated in section 3, the positivity of m?o constrains /i to lie in the range /i- < /i < /^h 



so as to satisfy Eq. (3.61). We give here the full expression for ii±: 



Xvj + 8(Ai + X,)v^ ± y^XjXvj + 16t;|((Ai + X,)vj + 4(A2 + X3)v!)) 



Note that due to the negative coefficient of /i in Eq.(3.61), ii± should always be real- 



of Eqs. (3.60 - 3.62), that is 



valued otherwise Eq.(3.61) is not satisfied and is tachyonic for all values of /x. As can 



be seen from Eq. (A.l), this requirement leads in principle to an extra constraint on top 



A(At;^ + 16t;2((Ai + X^)vl + 4(A2 + X3)vl)) > (A.2) 
However, we show here that this extra constraint is automatically satisfied due to the BFB 



constraints: since A > 0, cf. Eq. (6.1), it suffices to show that {Xvj + 16i;^((Ai + X4)v^ + 
4(A2 + Xs)v^)) > 0. Now using the first inequality of Eq. (6.3) one obtains 



Xvj + 16i;2((Ai + X,)vj + 4(A2 + Xs)v^) > (Xvj + IQv^- y^XiX^+X^jv} + 4(A2 + X,)v^)) 

> (VXvl - 8v/(A2 + Xs)v^)' (A.3) 



that proves our statement. 
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Appendix B. 

In this appendix we give the form of the Higgs potential in the field subspaces where only 
2 or only 3 fields are non-vanishing, dubbed respectively 2-field and 3-field directions. We 
identify exhaustively 10 different directions for each of these two classes and give their 
corresponding BFB conditions. 

The ten 2-field directions: 



T/{4) 




(B.l 


T/(4) 


= (A2 + A3) ir+l^ + (Ai + A4) ir+l^ + ^ 


(B.2 


T/(4) 


= (A2 + A3) + Ai i5++n0T + 7 


(B.3 


T/(4) 
2 l'dir.4 


= (A2 + ^) + (Ai + ^) + 7 


(B.4 


2 l^dir.S 


- (\^ + ^\ |A+|4 , ex , h\ |A+|2 U0|2 , A 1 ,0|4 


fB 5 


2 »^dir.6 


= (A2 + y) + 2 (A2 + A3) + (A2 + A3) \5+^\' 


(B.6 


2'/dir.7 


= (A2+A3)i5°r+Aii5Ti0+r+^i</'+r 


(B.7 


T/(4) 
2 >^dir.8 


= (A2 + A3) + (Ai + A4) i^Y i/r + ^ i/r 


(B.8 


T/(4) 
2^^dir.9 


= (A2 + A3) i^T + 2 A2 I'^^-'r + (A2 + A3) 


(B.9 


T/(4) 
2'^dir.lO 


= (A2 + A3) I^T + 2 (A2 + A3) |<5T + (A2 + y ) l^+r 


(B.IO 



direction 1 : A > (B.ll) 

directions 2 and 8 : A > 0, A2 + A3 > 0, Ai + A4 + VA(A2 + A3) > (B.12) 

direcitons 3 and 7 : A > 0, A2 + A3 > 0, Ai VA(A2 + A3) > (B.13) 

directions 4 and 5 : A > 0, A2 + ^ > 0, Ai + ^ + a / A(A2 + ^) > (B.14) 

A-, 

directions 6, 9 and 10 : A2 + A3 > 0, A2 + ^ > (B.15) 
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The ten 3-field directions: 



+ 2(A2 + A3)|5+H5++|^ + (A2 + A3)|5++|^ 

S- ,,0,2,. + ,4 . S'iS-'' 

50 (5-) 

3Kf^2 = (A2 + A3) \5r + 2(A2 + A3) |5T I'^'-l' + (A2 + ^) |5+ 



-A3 THTTZT^I^TI^I' - A3 ^4£r^ (B.16) 



2 

+ (Ai + A4) I^T I0T + (Ai + y ) l<^T + J (B.17) 



3Kf& = (A2 + A3) + 2(A2 + A3) |5T I^^P + (A2 + y) 15+'^ 



+x, \5r + (Ai + y) + ^ (b.is) 
3^r4 = (A2 + A3) |5T + 2 A2 |5T i^^+r + (A2 + A3) 

+ (Ai + A4) I^T 1/1' + Ai |5++r I0T + ^ 1/1' (B.19) 
sVaS., = (A2 + A3) \Sy + 2 A2 |<^++r + (A2 + A3) 

+ Ai|5ri</'+r + {Xi + X,)\S++\'\<p+\' + (B.20) 

3k;^6 = (A2 + A3) \5r + (Ai + A4) \5r i^t + ^ i<^t 

+Aii5Trr + ^i/ri0i' + ^i0i' (B.21) 

3Kf& = (A2 + ^) + 2(A2 + A3) \S^\' 15-^^1' + (A2 + A3) 

+ (Ai + y) |5+r I0T + Ai I0T + ^ I0T (B.22) 

3^r8 = (A2 + ^) + 2(A2 + A3) i<^-^-^r + (A2 + A3) 

+ (Ai + ^) |<5+r 101' + (Ai + A4) 101' + J 101' (B.23) 

3Kfi:!9 = (A2 + ^) 151^ + (Ai + ^) 151^ 10T + ^ i0r 

+ (Ai + ^) 101' + ^ I0T 101' + 7 101' (B-24) 



2 ' ' ' '"^ ' ' 2 '"^ ' '"^ ' 4 

A 

4 

A , ,n,r> , , I ,0 A 



3V^i- !io = (A2 + A3) + Ai |5+1^ I0T + 7 I0T 



+ (Ai + A4) ir+r 10+r + ^ I0T 10+1' + 2 I'^^l' (B-25) 
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The corresponding BFB conditions read: 



direction 1 : 2A2 + A3 > A A2 + A3 > A (^A^ < 4(A2 + A3)^ V A3 < (B.26) 
direction 2 : 

A > OA A2 + A3 > OA2A2 + A3 > A y/X{X2 + A3) + Ai + A4 > A (^(^2A(2A2 + A3) > (2Ai + A4)^ 

A (y2^Xs{X2 + A3) ((2Ai + A4)2 - 2A(2A2 + A3)) + 2A2 A4 > 2Ai A3 A A3 < o) V 

^(2A. + A3)((2A.+A.K2A.+3A.)-4A(A. + A3)) > „ A 2A, + A. < o) )) V 2A, + A. > o) (B.27) 



direction 3 



A > OA A2 + A3 > OA2A2 + A3 > A VA(A2 + A3) + Ai > A (^(^2A(2A2 + A3) > (2Ai + A4)^ A 
(2A2 + A3)(4A(A2 + A3)-4A2 + Al) . 

( {''^ + < ° ^ SaTTa^ < °) " 

((A2 + A3)(2A2 + A3 - 2) > A V2^Xs{X2 + A3) ((2Ai + X^)^ - 2A(2A2 + A3)) > 2A1A3 
+2A4(A2 + A3)) ) ) V 2Ai + A4 > 0) (B.28) 



direction 4 : 

A>0AA2 + A3>0A y/X{X2 + A3) + Ai + A4 > 

. (A2 + A3) (-AAi + A2(A2 - A3) + 2A1A2A4) 



A( ( ^ '->0 

A1A2 



a(^(^X2 > a A(A2 + A3) > A^ A Ai < 0) V (Ai > A A3(2A2 + A3) > A A2 < 0))) 
V(Ai > A A2 > 0) V (a(A2 + A3) > A? A A3(2A2 + A3) > 

A^A3(2A2 + A3) (a(A2 + A3) - Af) + A1A3 + A4(A2 + A3) > o)) (B.29) 



direction 5 



A>0AA2 + A3>0A v'A(A2 + A3) + Ai > A 

(A2 + A3) (aAI + XliXs - A2) + 2A1A3A4 + A2(A2 + A3)) 



< A 



A2(Ai + A4) 

(^(A3(2A2 + A3) > A Ai + A4 > A A2 < 0) V (^A2 > A A(A2 + A3) > (Ai + A4)^ 

AAi + A4 < 0))) V (A2 > A Ai + A4 > 0) V (a(A2 + A3) > (Ai + A4)^ A A3(2A2 + A3) > 

^J-X^{2X2 + A3) ((Ai + X^f - A(A2 + A3)) + Ai A3 > A2A4) ) (B.30) 



A 
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direction 6 : 

A>0AA2 + A3>0A ^/\{\2 + A3) + Ai > 

A (^Ai + A4 > V (^A(A2 + A3) > (Ai + A4)^ A A4 < 0) ) (B.31) 



direction 7 : 



A > A 2A2 + A3 > A A2 + A3 > A v/A(4A2 + 2A3) + 2Ai + A4 > A (^(^A(A2 + A3) > A? A 
((Ai(2A2 + 3A3) + 2A4(A2 + A3) > ^^(^^^+ A Ai < 0) V v^yA3(A2 + A3)(A2-A(A2 + A; 
+A4(A2 + A3) > 0)) VAi > 0) (B.32) 



direction 8 : 

A > A 2A2 + A3 > A A2 + A3 > A v/A(4A2 + 2A3) + 2Ai + A4 > A 

((a(A2 + A3) > (Ai + A4)' A ((v^yA3(A2 + A3) ((Ai + A4)2 - A(A2 + A3)) > A4(A2 + A3) A A3 < o) 
y (2X^X2 + 3A1A3 + A3A4 > "^^^^^^^^^^ A Ai + A4 < 0))) V Ai + A4 > 0) (B.33) 



direction 9 



A > A 2A2 + A3 > A V'A(4A^T2A^ + 2Ai + A4 > 

A (2X^2X2 + A3) > (2Ai + A4)^ V 2Ai + A4 > 0) (B.34) 



direction 10 : 

A>0AA2 + A3>0A VA(A2 + A3) + Ai + A4 > A (^Ai > V A(A2 + A3) > A^ V A4 > o) (B.35) 
We emphasize that all the above BFB conditions are contained in the general solution 



given by Eqs. (6.1 - 6.3). 



Appendix C. 

For completeness we give in subsections C.l and C.2 a partial list of the couplings in the 
DTHM that are relevant respectively to the discussion in section [8] and to the derivation 
of the results of section IH 

C.l Higgs-gauge boson couplings & triple Higgs couplings 



Shifting the neutral fields according to Eq. (2.26), and using the relations between the 
physical and non-physical state bases, 
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Z2 



(C.l) 



0± 
5± 



(C.2) 



whith TZa^Tlp and 7^/?/ as defined in Eqs.(2.12, 2.23), one extracts from tlie kinetic terms 



and tfie covariant derivatives, Eqs. ( |2.1 , 2.2, 2.3), tlie couplings involving Higgs bosons 



and gauge bosons, and from the potential, Eq. (2.4), the triple Higgs couplings 



We list below some of the resulting Feynman rules and provide also approximate ex- 
pressions in the limit of very small mixing between the triplet and doublet Higgs mul- 
tiplets, (i.e s„ = Oiy'^/vl), c^^^^p, = 1 + 0{v^/vj), = 2vt/vd + 0{vl/vl) and sp, = 
V2vt/v, + 0{v!/vj)). 



h'^ZZ = - 

H^ZZ = 
hPW+W 

h'A^Z = 
H'A^Z = 
h^h^H^ = 

h^W+H- 



g , 9 vt 
-t — mz{SaC,3 - '2caS(s)g^^ ^ 4z mzg^,u 

= igcwmz{caCp + SaSp)g^y ^ igmwgf^u 
- = -igmw{saC[^ - CaSfs)gf,u ~ 2igmw—g,M,y 

-T^iCaSp - 2Ci3Sa){Ph - Pa),i ~ — — {Ph- Pa),, 

2cw Cw Vd 

g \ g , 

7. [SaSp + 2CaCp)[pH - PAj^i ~ [PH - PAj^i 

Zcw Cw 
3 



(C.3 

(C.4 

(C.5 
(C.6 

(C.7 

(C.8 



CaVt) 



z ( V2/i + (3A - 5(Ai + K))vt ) + 0{vl) 



[CaSp> - V2SaCp'){ph-pH-)t, ~ +iA= — iPh-PH-),i 



\f2Vd 



■9 , 



A'^W+H- = ^-{V2cp'Cp + sp,sp){pA<^ 



Ph- 



g 

V2 



{Pao -PH-)t 



: gmz{cpSp'Sw 



SfjCjB 



-{l + slr))g,,^-V2g 



mzVt 



V2 " Vd 
H^+H'W^ = igCpi{pH++ -Ph-)^ ~ ig{PH++ - Ph-),, 
H^^W-W- = -iV2g\g^, 

H^^H^H^ = -z(2/is^, + C/3' (A4 Si3> Vd - V2\3 cp' Vt)) 
H++H—VV' = 8ievev'g,,u 
H++H-V = -2teviPH++-PH-),. 
H^H~VV' = 2ievev' g^^u 
H+H'V = -iev{pH+ -PH-)^i 
G+G-VV = 2ievev' g,,^ 
G^G'V = -iev{pG+ -PG-),i 



-g^iu 



(C.9 
(C.IO 

(C.ll 



(C.12 

(C.13 

(C.14 

(C.15 

(C.16 
(C.17 
(C.18 
(C.19 
(C.20 
(C.21 
(C.22 
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where in Eqs. (C.17 - C.22) we denote by V and V the 7 or Z gauge boson, with the 
couphngs satisfying = e and ez = ecot29w- We also adopted the convention that all 
momenta are incoming at each vertex. 



C.2 Quartic scalar couplings in the doublet-triplet basis 

Here we give the complete list of Feynman rules for the quartic scalar couplings in the 
unrotated basis which were used in section |5] to determine the unitarity constraints: 



6+ 6+ 6~ 6- = 
6+ 6- 6— 6++ 
5++ (5++ 6— 6- 
6+ 6+ 6— Z2 = 
6- 6- (5++ Z2 = 
6- (5+ Z2 Z2 = 
Z2 Z2 Z2 Z2 = 



-2z(2A2 + A3) , 

= -2t{X2 + X3) , 

- = -At{X2 + A3) 

-v^A3 , 
-22(A2 + A3) , 
-Qt{X2 + Xs) , 



Z2 Z2 
Zi 



6+ 6+ 6- 
6- 6- 6+ 
6+ Zi 0- 



0+ = -2(Ai) 
0+ = -ilX 
0+ = -tX 

- e = 

' e = ^V2Xs 

e 



Zi 0+ e 



A4 
-A4 

2^/2 



6++ 6— Z2 Z2 


= -2iX2 , 


6- 5+ e e = 


-2z(A2 + A3) 


6+ 6- Zi Zi = 


-|(2Ai + A4) , 


6— (5++ e e° 


= -2iA2 


S++ 6— Zi Zi 


= -^(Ai) , 


Z2 Z2 e e = 


-22(A2 + A3) 


Z2 Z2 Z\ Z\ = 


-z(Ai + A4) , 


Zx Zx e e = 


-^(Ai + A4) 


Zx Zx Zi Zx = 


-fa , 


0- 0+ ^0 ^0 = 


-z(Ai) 


S++ 6- Zi 0- 


_ A4 
2 ' 


^0 ^0 ^0 ^0 _ 


-6^(A2 + A3) 


(5+ 0- Zx Z2 = 


1X4 

2V2 ' 


6- 6++ 0- h -- 


_ *A4 
2 


5+ 6— 0+ Zx 


_ A4 
2 ' 


6+ (j)' Z2h = 


A4 
2V2 


6- 0+ Z2 Zx = 


1X4 

272 ' 


5+ 5— 0+ h -- 


_ 1X4 
2 


6- 5+ 0+ 0- = 


-|(2Ai + A4) , 


5" Z2 0+ /i = 


X4 
2V2 


6— 5++ 0+ 0- 


' = -t{Xx + A4) 


, S+ ^0 h(j)- = 


1X4 
2v^ 


6- h = - 


1X4 

2V2 ' 


Z2 Z2 h h = - 


-«(A4 + Ai) 


S- S+ hh = - 


i(2Ai + A4) , 


Zx Zx h h = - 




6— 6++ hh = 


-^(Ai) , 


(f)'^ (f)^ h h = 


^2 


hhi^ = -i 


(Ai + A4) , 


h h h h = -i^ 



One can read off from this list the /i-independent part of the triple scalar couplings, by 



making the substitutions Zx —ivd, h ^ or Z2 —ivt,C,^ — )■ Vf (cf. Eq. (2.26) ) in the 
appropriate vertices and modifying accordingly the symmetry factors for identical fields. 
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